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Abstract 

In  the  past  no  satisfactory  methods  have  been  available  for  determining 
the  reliability  of  one-shot  items  where  extremely  high  reliability  is 
required.  This  has  been  detrimental  to  the  development  of  one -shot, 
automatically  activated  batteries  for  missile  applications.  The  lack  of  a 
satisfactory  method,  one  which  will  permit  determination  and  stati stical 
analysis  of  the  failure  points  of  the  battery  with  respect  to  the  specified 
environmental  stresses,  has  been  compensated  by  reliability  programs  based 
on  designing  the  batteries  to  meet  environmental  stresses  with  large  safety 
factors. 

The  response  surface  determination  (RSD)  method  is  proposed  for  the 
determination  of  mean  failure  points  and  reliability  tolerance  limits  for  a 
battery  design  with  respect  to  operational  environmental  conditions,  thermal 
and  dynamic .  It  provides  for  the  analysis  of  the  battery  responses  as  a 
function  of  two  or  more  environmental  stresses  acting  simultaneously,  thus 
affording  information  on  the  effects  of  interactions  between  forces  on  the 
battery  performance .  Emphasis  is  placed  on  providing  a  maximum  of  reliabil¬ 
ity  prediction  data  with  small  test  sample  sizes. 
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RESPONSE  SURFACE  DETERMINATIONS  IN  ESTABLISHING 
THE  RELIABHITy  OF  ONE-SHOT  ITEMS 

INTRODUCTION 

This  report  deals  with  a  new  procedure  for  establishing  the  reliability 
of  one  -shot  items  under  simulated  operational  environne  nts  where  extremely 
high  reliability  is  required.  One-shot  items  are  herein  defined  as  components 
or  equipments  which  are  expended  in  use  and  which,  prior  to  use,  do  not  lend 
themselves  to  non -destructive  checkout  tests  from  which  their  probability 
of  successful  operation  or  reliability  can  be  inferred.  The  military  one- 
shot  items  of  perhaps  the  greatest  interest  to  USAELRDL  at  present  are  the 
automatically  activated  battery  power  supplies  which  have  been  developed 
for  guidance  and  control  functions  in  a  wide  range  of  missiles.  High 
reliability  standards  have  been  established  for  these  batteries  since 
failure  in  any  sense  will  destroy  the  effectiveness  of  the  missile. 
Unfortunately  no  effective  and  practical  method  has  been  available  to  date 
to  determine  the  probable  failure  rate  of  a  battery  design  under  simulated 
missile  environmental  conditions.  This  has  made  it  necessary  to  use 
compromise  measures  to  achieve  high  engineering  confidence  in  the  reliability 
of  the  battery  such  as  designing  it  to  meet  environmental  requirements  with 
high  Safety  factors  beyond  the  specified  requirements . 

The  environmental  safety  factor  concept1  has  resulted  in  a  commendable 
history  of  high  reliability  in  the  field  for  many  missile  batteries. 

However  it  has  been  challenged  from  many  sources  on  the  basis  of  unnecessarily 
increasing  the  weight,  size  and  cost  of  the  batteries.  The  concept  involves 
proving  the  battery  design  by  testing  relatively  small  samples  under  each 
critical  environmental  condition  where  the  g  force  for  shock,  vibration, 
acceleration,  spin,  etc.  is  held  at  some  arbitrarily  chosen  multiple  of  the 
required  dynamic  g  force,  normal ly  four  times.  Failure  of  the  battery  design 
in  any  manner  under  a  kX  environmental  force  results  in  redesign  until  the 
kX  capability  has  been  demonstrated.  Although  statistical  statements  cannot 
be  made  as  to  the  probability  of  failure  at  the  required  g  level,  the  concept 
does  provide  for  a  high  engineering  confidence  in  the  design. 

Other  approaches  to  battery  reliability  have  been  attempted  including 
statistical  ones  such  as  the  evaluation  of  the  variability  of  the  battery 
performance  parameters  [service  life,  activation  time,  maximum  voltage, 
minimum  voltage  under  heavy  load  pulses,  etc.]  under  bench  conditions  at  the 
required  temperatures  or  under  the  specified  dynamic  g  forces  or  even,  in  one 
program,  under  4X  g  forces.  The  statistical  approach  used  is  to  assume  a 
normal  distribution  for  the  battery  responses  and  then  to  evaluate  each 
battery  test  sample  for  the  mean  value  and  the  standard  deviation  for  each 
response.  The  number  of  standard  deviation  units  e.g.  between  the  mean 
service  life  and  the  required  service  life  must  exceed  a  specified 
reliability  standard,  or  K  factor^  which  is  related  to  the  sample  size,  the 
desired  maximum  failure  rate,  and  the  confidence  level  used  in  making  the 
reliability  statement .  Any  program  of  this  type  has  certain  basic  disadvan¬ 
tages.  No  real  insight  is  gained  into  the  effects  on  the  battery  performance 
of  the  thermal  and  dynamic  environmental  forces,  unless  this  is  sidestepped 
by  testing  at  the  h-X  levels.  Even  under  this  condition  the  testing  is  done 
with  one  environment  at  a  time,  thus  revealing  nothing  es  to  the  potentially 
significant  effects  of  interactions  between  the  environmental  variables . 
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It  was  recognized  for  some  time  that  the  environmental  variables  have 
a  major  effect  on  the  reliability  of  a  battery  and  that  some  test  method 
was  required  to  study  the  battery  performance  as  a  function  of  these 
environments  treated  as  continuous  variables.  In  fact  the  goal  was  to 
apply  the  same  type  of  statistics  as  discussed  before  to  demonstrate  that 
the  mean  failure  point  of  a  battery  with  respect  to  an  environment  would  be 
at  least  K  standard  deviation  units  above  the  required  stress  level.  A  first 
Step  in  this  direction  was  made  in  1961  by  the  development  of  a  step-stress 
technique3vhich  permitted  these  determinations.  The  method  had  certain 
basic  limitations  which  prevented  its  wide-scale  application.  It  assumed 
that  the  test  responses  were  normally  distributed.  It  provided  for  the  study 
of  only  one  environment  at  a  time.  Its  use  required  that  the  battery  be 
tested  at  stress  levels  high  enough  to  induce  failure.  This  proved  possible 
with  the  thermal  environments,  high  and  low  temperature,  but  it  was  rapidly 
shown  that  the  method  was  not  generally  adaptable  to  dynamic  environments 
since  practical  test  equipment  in  many  if  not  most  cases  could  not  induce 
failure. 

The  RSD  method  has  been  developed  at  USAELRDL  during  the  past  two  years 
to  overcame  these  disadvantages.  Although  it  has  not  been  applied  as  yet 
in  an  actual  reliability  program,  it  is  believed  capable  of  providing  the 
desired  environmental  capability  data  for  a  given  battery  design.  Determin¬ 
ing  the  mean  minus  KB  reliability  tolerance  limit  should  be  possible  without 
testing  to  failure.  Within  the  limits  of  practical  test  equipment,  the 
method  should  provide  for  the  testing  under  two  or  more  different  environ¬ 
mental  stresses  acting  simultaneously, •  thus  affording  some  of  the  desired 
insight  into  potentially  destructive  interactions  between  environmental 
stresses.  Instead  of  merely  assuming  normal  distribution  of  responses, 
the  method  provides  a  test  which,  though  not  conclusive,  at  least  provides 
assurance  in  applying  normal  distribution  statistics  in  making  reliability 
predictions. 

DISCUSSION 

The  Response  Surface  Determination  (BSD)  Method 

The  basic  techniques  of  the  RSD  method  are  based  on  conventional  design 
of  experiments  procedures  which  are  thoroughly  discussed  in  the  literature ,4 
These  procedures  provide  for  the  study  of  the  response  of  an  item  as  a 
function  of  two  or  more  continuous  variables  acting  simultaneously.  They 
are  normally  used  to  provide  information  as  to  the  settings  of  the  variables 
to  yield  an  optimum  condition  for  the  response. 

Translated  into  battery  testing  terms,  the  continuous  variables  of 
interest  are  the  thermal  and  dynamic  environments  to  which  the  battery  is 
exposed  in  the  missile  operation.  The  responses  to  be  evaluated  are  the 
battery's  critical  performance  parameters:  activation  time,  service  life 
or  capacity  to  minimum  voltage,  and  maximum  and  minimum  voltages  under 
continuous  or  pulsating  loads.  Each  of  these  responses  are  evaluated 
separately  from  the  same  test  data  and  individual  reliability  predictions 
are  made  for  each  response. 


2 


The  conventional  design  of  experiments  procedures  give  an  indication  of 
the  mean  or  average  response, hut  they  do  not  provide  for  tolerance  intervals 
based  on  reliability  standards.  The  literature  only  covers  the  simple  case 
of  tolerance  intervals  around  a  linear  regression  line  involving  one 
variable.  There  is  no  practical  technique  given  for  curved  functions 
involving  two  or  more  variables.  Straightforward  techniques  were  therefore 
derived  for  this  and  these,  together  with  the  design  of  experiments  procedures, 
constitute  the  PSD  method.  This  report  is  concerned  entirely  with  the  study 
of  the  battery  response  as  a  function  of  only  two  continuous  variables  since 
present  testing  equipment  for  the  batteries  permits  only  two  environments, 
one  thermal  and  one  dynamic,  to  be  applied  simultaneously.  For  one-shot 
items  capable  of  testing  under  more  than  two  environments  simultaneously, 
or  with  the  advent  of  more  sophisticated  test  equipment,  the  techniques  as 
given  will  have  to  be  expanded.  This  is  entirely  possible  and  consistent 
with  conventional  design  of  experiments  procedures,  though  beyond  the  scope 
of  this  report. 

The  PSD  theory  will  be  developed  through  a  series  of  illustrations . 

The  mechanics  of  the  operations  will  be  developed  later  through  a  series  of 
examples.  Figure  1  shows  a  three  dimensional  space  generated  by  the  Y  or 
response  axis  (e.g.,  battery  capacity) ,  the  X-j_  axis  (e.g.,  high  temperature 
as  a  continuous  variable )  and  the  Xg  axis  (e.g.,  vibration  g  force,  also 
as  a  continuous  variable ) .  Point  A,  the  intersection  of  the  axes, 
represents  the  combined  requirement  point,  e.g.,  10  seconds  capacity  at  a 
10  g  vibration  force  at  165  °F.  The  objective  is  to  determine  that  the 
battery  has  the  desired  degree  of  reliability  with  respect  to  capacity  at 
point  A  and  also  to  determine  how  far  along  the  X  axes  the  given  reliability 
standards  Will  be  met,  and  also  for  any  given  combination  of  X^  and  Xg  above 
the  requirement  levels. 

The  first  step  before  exploring  the  three  dimensional  space  is  to 
determine  if  the  battery  design  will  meet  the  reliability  standards  at 
point  A  (Figure  2).  A  sufficient  sample  is  tested  at  point  A,  10  g  and 
165  °F,  and  the  data  is  analyzed  to  determine  the  mean  capacity  Y  and  the 
standard  deviation,  or  measure  of  variability  of  individual  points  around 
the  mean.  A  test  is  made  to  see  if  it  is  reasonable  to  assume  that  the 
points  are  normally  distributed.  (The  techniques  for  all  of  these  procedures 
will  be  ejq)lained  in  detail  in  the  later  experimental  section) .  If  the 
assumption  of  a  normal  distribution  can  be  made,  a  normal  curve  can  be 
superimposed  on  the  Y  axis,  laid  off  in  standard  deviation  units,  and  the 
Y  -  KS  capacity  value  can  be  determined,  where  K  represents  the  reliability 
standards ?  Any  Y  -  K5  value  over  10  seconds  indicates  adequate  reliability 
at  point  A. 

As  shown  in  Figure  3>  the  XqXg  test  design  is  then  laid  out.  One  or  two 
batteries  are  first  run  along  the  X^  and  Xg  axes  to  obtain  an  approximate 
idea  of  the  mean  failure  temperature  and  g  force  levels .  As  will  often 
prove  the  case,  the  upper  limit  of  the  vibration  equipment  capability  will 
not  induce  a  failure .  In  this  case  this  limit,  hereafter  referred  to  as 
the  test  equipment  capability  or  TEC,  will  form  the  Xg  upper  boundary.  Some 
temperature  B  will  form  the  boundary  along  the  Xj  axis.  Points  C  and  D  are 
established  at  the  midpoints  of  the  A-B  and  A -TEC  intervals .  This  is  a  3^ 
factorial  design  with  nine  test  points,  eight  of  which  remain  to  be  investi¬ 
gated.  Any  appropriate  scale  can  be  used  along  each  axis  as  long  as  equal 
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increments  are  held.  In  the  case  of  vibration,  assuming  that  the  TEC  is 
ho  g,  the  midpoint  D  can  be  25  g  -with  a  linear  scale  or  26  g  with  a 
logarithmic  Scale,  equal  increments  in  both  cases.  Assuming  B  to  be  205 °P 
C  becomes  185®F. 

Two  tests  are  then  conducted  at  each  of  the  eight  new  test  points 
giving  capacities,  Y  responses,  as  shown  on  Figure  k.  The  test  data  points, 
as  shown,  may  be  considered  typical  of  what  might  be  expected  if  the  battery 
response  is  influenced  by  the  environments.  Looking  at  the  points,  at 
10,  26  or  1*6  g  the  dropoff  or  regression  of  Y  with  respect  to  definitely 
does  not  appear  to  be  linear  but  appears  to  suggest  a  mathematical  relation¬ 
ship  involving  second  order,  quadratic,  terms.  Similarly  with  the  X^ 
points.  The  depression  in  the  responses  with  increasing  levels  of  XT  gets 
more  pronounced  at  higher  levels  of  Xg.  The  Same  applies  in  the  other 
direction.  This  is  the  interaction  effect  which  means  that  the  effect  on 
the  response  of  one  variable  cannot  be  stated  properly  without  referring  to 
the  level  of  the  other  variable.  Looking  at  the  points  as  a  whole  it  can 
be  seen  that  ajcurved  response  surface  can  be  fitted  through  the  points 
starting  with  Y^  and  dropping  off  with  increasing  levels  of  X^  and  Xg. 

This  surface  is  shown  on  Figure  5. 

The  response  surface  does  not  fit  uniformly  through  the  test  points 
but  instead  has  a  unique  property.  It  is  technically  known  as  a  least 
squares  regression  Surface,  which  means  that  if  the  vertical  deviations  of 
the  18  test  points  (doubling  Y^)  from  the  surface  are  squared,  then  the  sum 
of  these  squares  will  be  at  a  minimum  value  for  this  one  Surface.  It  is 
based  on  a  two  variable,  second  order  mathematical  model  of  the  form: 

Y  ■  bg  +  b^X^  +  bgXg  +  b-[  i  X-j  ^  +  bppXg^  +  b^pX^Xg 

where  the  coefficients  are  derived  from  analysis  of  the  l8  test  points. 

The  X^Xg  term  represents  the  interaction.  The  response  surface  is  an 
estimate  of  the  mean  value  of  Y  for  any  settings  of  X.  and  Xg  throughout 
the  test  area.  Therefore  the  intersection  of  the  surface  with  the  16  second 
Y  plane  is  the  mean  failure  contour  which  represents  an  important  estimate 
of  the  population  from  which  the  battery  sample  was  drawn.  The  estimate 
is  that  if  additional  samples  were  taken,  half  would  pass  and  half  would 
fail  at  any  X^  and  X2  settings  along  the  contour.  But  this  is  of  little 
value  in  predictingjreliability  of  the  battery  design.  What  is  required 
is  another  surface  Y  -  KS  below  the  response  surface,  the  intersection  of 
which  with  the  10  second  plane  will  be  the  reliability  boundary  representing 
the  highest  values  of  settings  of  X^  and  Xg  at  which  the  reliability 
standards  can  be  met.  A  parallel  surface  drawn  through  (Y  -  KS)^  would  not 
suffice  since  there  is  no  reason  to  assume  that  this  estimate  of  variance 
applies  throughout  the  whole  test  region.  Likewise  it  isn't  known  that  the 
variance  is  uniform  throughout  and  that  the  responses  are  normally 
distributed  around  the  response  surface  throughout  the  whole  test  region. 

An  indication  of  a  normal  distribution  had  been  demonstrated  at  the 
16  g,  165 °F  point.  A  similar  accumulation  of  data  is  not  available  elsewhere 
in  the  test  region  to  make  an  additional  test  for  normality.  Therefore  a 
basic  assumption  is  made  that  if  normality  can  be  assumed  at  the  combined 
requirement  point,  it  can  be  assumed  throughout  the  entire  test  region.  It  ia 
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recognized  that  if  departures  from  normality  occur,  they  will  influence  the 
position  of  the  reliability  boundary,  possibly  weakening  the  value  of  these 
predictions.  However  they  will  not  Weaken  the  value  of  the  prediction  at 
the  combined  requirement  point,  the  most  important  of  the  reliability 
predictions.  The  prediction  as  to  how  high  we  can  go  in  temperature  and 
still  achieve  reliable  operation  may  be  influenced  slightly,  but  this  is  not 
as  critical  a  prediction.  The  assumption  of  uniform  variance  can  be  tested 
for  with  the  available  data.  This  is  done  by  determining  the  variance 
throughout  the  test  area  (except  for  the  combined  requirement  point)  and 
comparing  it  to  the  variance  previously  determined  at  the  combined  requirement 
poinx.  If  these  two  estimates  of  the  variance  cannot  be  proved  to  be 
significantly  different,  then  uniform  variance  can  be  assumed  by  inference, 
at  not  too  great  a  risk.  If  uniform  variance  can  be  assumed  in  this  manner, 
the  two  estimates  of  variance  can  then  be  pooled  to  give  a  common  estimate, 
based  on  all  of  the  test  data,  of  the  variance  and  from  this  of  the  standard 
error  throughout  the  entire  test  region.  This  latter  estimate,  multiplied 
by  K,  is  then  used  to  establish  a  Y  -  KS  equation  for  a  reliability  boundary 
surface  parallel  to  the  response  surface  as  shown  in  Figure  6. 

The_reliability  boundary  surface  is  drawn  parallel  to  the  response 
surface,  Y  -  KS  below  it.  The  intersection  with  the  10  second  plane  is  the 
required  reliability  boundary.  The  basic  estimate  now  made  from  the  Sample 
is  that  the  reliability  standards  will  be  met  at  any  combination  of  and 
X2  to  the  left  of  the  reliability  boundary  including,  of  course,  the  combined 
requirement  point  of  l65°F  and  10  g.  No  statement  should  be  made,  by 
extrapolation  of  the  reliability  boundary,  of  how  many  g  it  would  be  possible 
to  go  to  at  165 °F  and  still  meet  the  reliability  standards.  The  only  valid 
statement  of  this  is,  in  excess  of  ^0  g.  The  interrelation  between  the  X 
and  Y  axes  should  be  noted.  Predictions  can  be  made  for  different  settings 
of  Y  as- well.  An  11  second  plane  could  be  added,  for  example,  and  the 
intersection  of  this  plane  with  the  reliability  boundary  surface  would 
indicate  the  maximum  settings  of  X^  and  Xg  at  which  11  seconds  of  operation 
could  be  achieved  in  accordance  with  the  same  reliability  standards.  A 
limitless  amount  of  reliability  prediction  data  could  actually  be  obtained 
once  the  reliability  boundary  surface  has  been  determined,  a  further 
advantage  of  this  approach  to  reliability. 

The  problem  arises  as  to  how  positively  the  reliability  predictions  can 
be  made.  If  there  were  no  doubt  at  all  that  the  responses  were  normally 
distributed  about  the  response  surface  throughout  the  entire  test  region, 
and  that  the  variance  were  completely  uniform  throughout  the  region,  then 
there  would  be  no  hesitancy  in  stating  that  we  are,  for  example, 
confident  in  predicting  that  no  more  than  one  battery  in  a  thousand  will 
fail  when  tested  at  193°F  while  undergoing  a  10  g  vibration.  However  the 
sample  on  which  this  estimate  was  based  was  also  used  to  provide  assurance 
in  assuming  normal  distribution  and  uniform  variance.  It  follows  logically 
that  the  size  of  the  sample  must  bear  heavily  on  the  faith  in  the  prediction. 
Further  comments  will  be  made  on  this  subject  in  the  Conclusions  section  of 
this  report  once  the  mechanics  of  the  BSD  method  have  been  explained  and 
the  area  of  sample  size  has  been  explored. 

Tbmp-H  mental.  Procedure 


The  preceding  section  has  attempted  to  explain  the  basic  principles  of 
the  BSD  method  and  general  considerations  governing  its  use  in  a  practical 

10 


RELIABILITY  BOUNDARY 
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test  program  for  determine  the  operational  reliability  of  one-shot  items 
Such  as  battery  power  supplies  for  missiles.  However,  in  this  area,  as  in 
most  in  the  field  of  applied  statistics,  no  intuitive  understanding  is 
possible  for  the  practicing  engineer  without  detailed  examples  to  study. 
Extensive  examples  have  been  prepared.  These  will  be  introduced  by  taking 
what  is  believed  to  be  representative  test  data  which  will  be  treated, 
following  the  RSD  method,  in  a  step  by  step  fashion  through  to  the  final 
reliability  conclusions.  Pull  details  of  all  calculations  are  given  in 
sequence  in  Appendix  I. 

!Exgeriaent_^_ 

Among  its  other  requirements,  a  battery  is  required  to  deliver  a 
minimum  of  10  seconds  service  life  within  specified  voltage  tolerances  at  a 
specified  maximum  temperature  of  165  °F,  while  being  subjected  to  a 
sinusoidal  vibrational  force  of  lOg  throughout  the  specified  frequency 
range.  It  is  desired  to  demonstrate  that  the  battery  design  has  a 
reliability,  or  probability  of  successful  operation,  of  a  least  99*95 » 
expressed  with  a  confidence  level  of  95$*  i*e. ,  admitting  a  chance  of  only 
one  in  20  of  stating  that  the  battery  has  this  reliability  when,  in  fact, 
it  does  not.  It  is  also  desired  to  know  the  maximum  temperature  at  which 
the  battery  will  operate  reliably  when  being  vibrated  at  10  g,  and  also  the 
maximum  g  force  to  which  the  battery  may  be  subjected  while  still  main¬ 
taining  95$  confidence  in  predicting  that  only  one  battery  out  of  a 
•thousand  will  give  less  than  10  seconds  of  operation  when  discharged  at 
165  F.  Out  of  the  total  quantity  of  batteries  available  in  the  qualification 
test  lot,  it  is  desired  to  demonstrate  the  reliability  with  as  small  a  test 
sample  as  possible  to  assure  sufficient  batteries  being  available  for  the 
study  of  other  environmental  variables  which  may  prove  to  be  more  detrimental 
than  the  present  high  temperature -vibration  combination.  As  a  starting 
point,  a  small  test  sample  of  15  Is  drawn  at  random  from  the  lot. 

Step  1.  Seven  batteries  are  selected  for  discharge  test  at  the 
composite  requirement  point  of  165  F  (X-^)  and  10  g  (Xg),  leaving  eight  for 
later  exploration  of  the  other  eight  points  of  the  X^,  Xg  space  of  the  3^ 
test  design.  The  batteries  are  tested  one  at  a  time  by  stabilizing  at  165  F 
and  then  quickly  removing  them  from  the  ambient,  mounting  them  on  the 
vibration  equipment,  activating  and  discharging  them  before  any  cooling 
effect  takes  place  (or  with  appropriate  temperature  monitoring  techniques 
to  assure  accurate  temperature  control).  The  batteries  yield  the  following 
service  times  in  seconds  (Y  responses)  to  minimm  voltage: 

15-74 
15.63 
ik.  in 
16.61 
15*47 
14.51 
18.79 

Stejj^2.  The  results  are  analyzed  (appendix  Ia)  for  the  mean  capacity, 
15.8b  seconds,  and  the  standard  deviation,  1.488  seconds.  The  reliability 
defining  K  factor  for  a  sample  of  size  H*7>  99*9$  reliability,  and  95$ 
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confidence  level  is  6.06f.  Therefore,  the  critical  Y  -  KS  value  is  15*88  - 
6.06l  x  1.488  ■  6.86  seconds.  This  is  well  below  the  10  second  requirement. 
This  indicates  that  either  the  battery  design  is  basically  unreliable  or 
that  perhaps  the  sample  size  was  too  small  to  allow  the  design,  if  good,  to 
reveal  itself  as  such.  This  latter  view  may  be  considered  valid  Since  the 
primary  objective  is  to  concede  the  unreliability,  and  to  then  search  for 
reliability  tolerance  limits  below  the  requirement  levels,  only  when 
convinced  that  the  reliability  standards  cannot  be  met.  Evidence  for  Such 
a  decision  may  be  based  on  trends  in  the  Y  *  KS  level  with  additional 
testing,  up  to  a  practical  limit  of,  for  example,  20  or  25  units.  -^Phis 
approach  will  now  be  explored  for  the  present  example  by  testing  additional 
randomly  selected  units  one  at  a  time  and  observing  the  trend  in  Y  -  KS  with 
each  sample  as  tested.  The  following  results  are  obtained: 


N 

Y 

Y 

K 

S 

KS 

Y  -  KS 

“7 - 

_ 15.6b _ 

6.061 

1.466 

9.02 

6.66 

8 

14.40 

15.70 

5.686 

1.474 

8.38 

7.32 

9 

17.66 

15.91 

5.4l4 

1.526 

8.26 

7.65 

10 

16.43 

15.97 

5.203 

1.442 

7.50 

8.4? 

11 

16.38 

16.00 

5.036 

1.380 

6.95 

9.05 

12 

16.724- 

16.06 

4.900 

1.333 

6.53 

9.53 

13 

15.61 

16.03 

4.787 

1.282 

6.14 

9.89 

14 

14.33 

15.91 

4.690 

1.313 

6.16 

9.75 

15 

17.52 

16.02 

4.607 

1.332 

6.14 

9.88 

16 

16. 04 

16.02 

4.53*4 

1.287 

5.84 

10.16 

17 

17.73 

16.12 

4.471 

1.313 

5.87 

10.25 

Graphically,  the  trend  is  shown  in  Figure  7.  In  this  example,  the  trend 
is  quite  clearly  upward  due  to  a  slight  increase  in  Y  and  a  definite  decrease 
in  S.  The  testing  is  stopped  according  to  an  arbitrary  rule  that  two 
successive  Y  -  KS  values  above  the  capacity  requirement  will  establish  the 
design  reliability  with  an  upward  trend  such  as  this.  The  reduction  in_K 
with  increasing  sample  size  is,  of  course,  a  factor  in  the  increase  of  Y  -  KS, 
but  not  a  significant  one  if  the  design_actually  had  no  chance  of  establishing 
itself.  To  indicate  this,  assume  that  Y  and  S  remained  at  the  N  =  7  levels 
of  15.88  and  1,2*88.  The  increase  in  Y  -  KS  with  increasing  N  is  as  shown 
by  the  lower  line  which  eventually  reaches  10  seconds  at  N  =  35.  In  practice, 
the  existence  of  an  upward  trend,  when  testing  up  to  a  practical  limit  of 
20  or  25  units,  should  be  judged  by  an  apparently  significant  difference  in 
the  two  lines,  as  shown  by  the  shaded  area.  With  an  inherently  unreliable 
design,  testing  could  be  stopped  presumably  after  a  total  N  of  10  or  15 . 

In  any  case,  if  the  required  Y  -  IS  value  is  not  reached  by  IT  -  20  or  25, 
testing  is  stopped  and  the  response  surface  test  design  is  then  centered  on 
the  Xj_,  Xg  requirement  point. 

StepS.  Having  judged  the  battery  design  to  have  adequate  reliability 
at  the  &L,  Ikg  requirement  point,  the  next  step  is  to  make  a  determination  of 
the  approximate  normality  of  the  distribution  of  the  individual  test  result 
points  about  the  mean  value  of  16.12  for  N  =  17.  This  is  done  by  plotting 
the  cumulative  distribution  of  the  sample  responses  on  arithmetical 
probability  paper  and  observing  if  an  approximate  straight  line  is  formed 
by  the  points,  particularly  in  the  range  of  P  «  10  to  P  »  90*  Calculations 
for  the  plotted  points  are  given  in  Appendix  IB.  The  plot  is  given  in 
Figure  8. 
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Trend  in  Y  —  KS  with  Increasing  N 


Cumulative  Freauency  Plot 


The  points  are  noted  to  fall  along  a  straight  line  very  well.  A 
straight  line  is  the  plot  one  obtains  when  plotting  the  cumulative 
distribution  of  a  large  or  infinite  Sample  from  a  known  normally  distributed 
population.  Therefore,  it  is  highly  probable  that  the  test  sample  was 
drawn  from  a  normally  distributed  population.  If  the  population  were  not 
normal i y  distributed,  in  which  case  it  would  plot  out  as  some  curved  line, 
then  the  probability  is  that  the  test  sample  plot  would  deviate  from  the 
straight  line. 

This  test  is  regarded  as  an  empirical  indication  of  the  suitability  of 
assuming  a  normal  distribution  of  individual  Y  responses  about  the  mean 
value  at  the  combined  X^,  Xg  requirement  point.  Having  satisfied  the 
assumption  at  the  most  important  point  of  the  X^,  Xp  space,  normality  is 
now  assumed  throughout  the  entire  test  range  of  the  Xg  space,  recognizing, 
of  course,  that  any  departures  from  normality  may  have  an  "influence  on  the 
eventual  estimates  of  the  battery  reliability  capabilities  along  the  X^  and 
Xg  axes,  but  not  on  the  estimates  in  the  region  close  to  the  X^,  X^  require¬ 
ment  point. 

Step  to  Having  demonstrated  reasons  for  assuming  the  reliability  and 
the  normal 1 ty  at  the  combined  X, ,  Xg  requirement  point,  and  having  obtained 
an  estimate,  16.12,  of  the  population  mean  and  an  estimate,  1.725,  of  the 
population  variance  (S2)  at  that  point,  the  test  design  to  explore  the  X^, 

Xg  space  will  be  established.  This  will  be  a  32  factorial  design  with  its 
lower  boundaries  established  by  the  X.  and  Xg  requirement  levels.  The  upper 
boundaries  Will  be  established,  as  explained  in  the  previous  section, 
preferably  at  an  estimate  of  the  mean  failure  point  but,  where  this  lies 
above  the  TEC  level  for  a  given  environment,  then  at  the  TEC  level.  In  this 
example,  let  it  be  assumed  that  the  mean  failure  point  at  Xg  *  10  g  has  been 
estimated  at  205 *F.  Let  it  also  be  assumed  that  the  mean  failure  point  at 
165 ®F  is  in  excess  of  to  g,  the  vibration  TEC.  Therefore,  the  test  design 
is  established  with  X^  limits  of  165  °F  and  205 °F  and  Xg  limits  of  10  g  and 
to  g.  The  intermediate  X.  value  is  l85°F.  The  intermediate  Xg  value, 
ins+ead  of  being  (to  +  10/2  m  25  g,  is  set  at  20  g  to  illustrate  the  technique 
of  transforming  a  variable,  in  this  case,  on  a  logarithmic  basis  in  which 
case  10  to  20  and  20  to  to  g  represent  equal  increments  along  Xg.  Coded 
values  are  assigned  to  the  X^  and  Xg  values  as  -1,  0  and  1  levels.  These 
coded  values  are  of  the  utmost  importance  for  keeping  calculations  as  simple 
as  possible .  All  equations  and  expressions  from  now  on  in  this  example  will 
be  based  on  the  coded  values.  Thus,  the  combined  X^,  Xg  requirement  points 
will  be  referred  to  as  (-1,  -l),  the  185 °F,  to  g  point  as  (0,1),  etc.  Based 
on  the  above,  the  test  design  is  as  shown  in  Figure  9* 

Stegjj.  Tests  are  now  run  at  each  of  the  eight  points  of  the  design 
otherTEan  ( -1,  -l).  The  results  are  as  follows: 
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The  responses,  shown  on  the  test  design,  are  as  follows.  With  the  16.12 
mean  value  introduced  for  the  ( -1,  -1)  point: 


*2 


The  responses  indicate  potentially  significant  effects  on  the  response 
with  increasing  stress  levels  of  X^  and  Xg.  The  coefficients  are  obtained 
(Appendix  IC)  for  the  two  dimensional,  second  order  mathematical  model  for 
the  response  surface: 

y  *  13.52  *  3.67X3^  - 1.673^  -  0.76X3*  -  0.20X2*  -  1.13X3X2 

An  analysis  of  variance  (Appendix  ID)  is  now  conducted  to  test  the 
equation.  Since  the  tests  were  not  replicated  (repeated)  at  the  eight  test 
points,  no  measure  of  the  experimental  error  can  be  made  and  therefore  no 
measure  of  the  "lack  of  fit"  of  the  model  to  the  data  can  be  made.  The  only 
important  test  of  the  data  that  can  be  made  is  the  determination  of  whether 
or  not  the  second  order  coefficients,  b^,  b—  end  b^,  are  significant. 

It  is  seen  from  the  analysis  of  variance  that  it  cannot  be  concluded,  that 
the  quadratic  components  are  significantly  different  from  zero.  Therefore 
it  may  be  concluded  that  a  linear  model  of  the  form 

*  -  b0  ♦  \\  +  V2 

would  suffice.  However  in  this  case  such  a  model  is  not  fitted  because 
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there  is  reason  to  doubt,  for  theoretical  reasons,  that  the  quadratic  terms 
are  not  significant.  Therefore  the  test  design  is  replicated  at  the  eight 
points  (except  for  (-1,  -l))  to  check  this  as  well  as  to  check  the 
experimental  error  and  the  "lack-of -fit"  term.  The  additional  test  results, 
added  to  those  previously  determined,  are: 


The  average  responses,  shown  on  the  test  design,  are  as  follows  With 
the  16.12  mean  value  again  introduced  for  the  ( -1,  -l)  point. 


The  coefficients  are  obtained  (Appendix  IE)  for  the  two  dimensional,  second 
order  mathematical  model  based  on  the  replicated  test  data: 

y  „  ik. 33  -  3.7OX1  -  1.70X3  -  1. 65X^2  -  0.60X22  -  0.92X3X2 

An  analysis  of  variance  (Appendix  IF)  is  now  conducted  to  test  for  the 
significance  of  the  quadratic  components,  b-, ,  bgg  and  bjp.  Due  to  the 
existence  of  the  extra  test  data,  it  is  possible  to  divide  the  residual 
sum  of  squares  into  two  parts,  that  due  to  experimental  error  and  whatever 
is  left  which  is  called  the  lack  of  fit  term,  a  measure  of  the  inability 
of  the  mathematical  model  to  fit  the  test  data. 

It  is  found  that  the  quadratic  terms  are  significant  and  that  the  lack- 
of  -fit-  term  is  not  significant.  Therefore  the  mathematical  model  is 
accepted  as  an  adequate  interpretation  of  the  experimental  data. 
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_Step6^  it  is  now  necessary  to  study  the  variance  throughout  the  test 
Space  and  to  compare  it  to  the  variance  at  the  (-1,  -l)  requirement  point. 

The  variance  is  determined  (Appendix  IG)  by  summing  the  squares  of  the 
differences  between  each  of  the  16  test  results  and  the  values  predicted  at 
the  test  points  from  the  equation.  The  variance  is  found  to  be  1.00893, 
which  is  to  be  compared  to  the  variance  estimate  of  1.72456  found  for  the 
17  results  at  (-1,  *1).  The  test  region  variance,  1.00893*  ia  based  on 
N  -  5  *  11  degrees  of  freedom.  The  (-1,  -l)  variance  was  based  on  N  -  1  ■ 

17  -1  *  16  d.f .  The  Fn  16  ratio  -  1.00893/1.72456  *  0.59.  This  ratio 
falls  within  the  limits  if  f  «  .  of  2.94  and  of  1 

vll,  lo  d.f.;  '■“■!■»  i'i  .1  l  m  . 

F . 95  f  ( lb,  11  d.f . ) 

of  0.302.  These  are  the  95$  level  of  significance  limits  of  an  F  test^for 
determining  acceptance  or  rejection  of  the  hypothesis  that  there  is  no 
Significant  difference  between  the  two  estimates  of  variance.  In  thi6  case 
the  hypothesis  is  accepted  that  there  is  no  significant  difference  between 
the  variance  at  (-1,  *1)  and  the  variance  throughout  the  remainder  of  the 
test  area.  It  may  therefore  be  concluded  by  inference  that  there  is  a 
uniform  variance  throughout  the  entire  test  region.  The  two  estimates  of 
variance  are  then  pooled  to  give  an  estimate  of  S,  the  standard  error 
throughout  the  entire  test  region.  This  value,  from  Appendix  m,  is  1.197 
seconds. 

_Steg7.  The  total  sample  size  used  in  this  exploration  of  the  X-jXg 
space  was  33,  17  at  (-1,  -1)  and  two  at  each  of  the  other  eight  points.  Of 
the  starting  33  degrees  of  freedom,  six  were  used  in  establishing  the 
coefficients  of  the  mathematical  model,  and  the  remaining  27  for  establishing 
the  pooled  standard  error.  Twenty-seven  is  also  the  effective  sample  size 
in  determining  the  K  factor,  which  for  N  °  27,  P  *  0.999  and  Confidence  = 

0.95  is  4.090.  The  KB  factor  is  therefore  4.090  x  1.197  ■  4.90  seconds. 

This  factor  is  used  in  establishing  the  Y  -  KS  reliability  boundary  value 
for  any  point  in  the  X^,  Xp  test  space.  Actual  equations  may  now  he 
established  for  finding  values  at  any  settings  in  the  three  dimensional 
Y,  Xj_,  Xg  space.  Two  of  the  most  convenient  equations  are: 

1.  The  equation  for  mean  failure  points  of  Y  with  respect  to 

Xj_  and  Xq’. 

Y  «  14.33  -  3.70X1  -  1.70X3  -  1.65X-L1 2  -  O.6OX32  -  0. 92X3X3  = 

10  (Y  requirement)  or  3.70Xi  +  1.70X3  +  1.65X12  +  0.60Xg2  + 
O.92X3X3  -  4.33 

2.  The  equation  for  reliability  boundary  points  of  Y  with  respect 
to  Xi  euad  Xg: 

7  -  14.33  -  3.70X3  -  1.65X12  -  O.6OX32  -  0.92XLX2  -  10  +  4.90(KS)« 

14.90  or  3.70Xi  +  1.70Xg  +  I.65X32  +  0.60X^2  +  0.92XjXg  = 

-0.57. 

Substituting  appropriate  values  of  Xi  or  Xg  (coded  values,  not  actual  ones) 
and  solving  the  resulting  quadratic  equations  (Appendix  XX  and  U) 
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values  from  which  the  mean  failure  contour  (l)  and  the  reliability  boundary 
contour  (2)  are  plotted  in  Figure  9.  (The  significance  of  curves  (A)  and  (B) 
will  he  explained  later). 

Analj^is^of^FxDeriment^ 

The  experiment  is  substantially  completed  with  the  plotting  of  the  mean 
failure  contour  and  the  more  significant  reliability  boundary.  Much  more 
could  of  course  be  done  if  necessary  or  desired.  For  example,  if  the 
reliability  boundary  had  passed  below  the  (-1,  *1)  requirement  point 
(indicating  a  probability  of  more  than  one  unit  per  thousand  failing  to  give 
10  seconds  of  service  at  (-1,  *1),  expressed  with  a  95#  confidence  level) 
then  instead  of  launching  a  redesign  of  the  battery,  the  missile  user  may 
decide  that  he  will  be  satisfied  if  he  can  reliably  get  nine  seconds  of 
operation.  An  equation  (3*70%  +  l,70Xg  +  1.6$%^  +  0.60Xg^  +  0.92X2%  * 
-1.57)  will  then  permit  determination  of  the  nine  second  reliability 
boundary  (intersection  of  the  surface,  which  is  Y  -  KB  below  the  mean 
response  surface,  and  the  nine  second  plane).  For  a  case  of  adequate 
reliability  such  as  the  present  example,  approximations  of  the  X^  and  % 
reliability  boundaries  can  be  made  for  higher  requirement  levels  of  Y,  11 
or  perhaps  even  12  seconds.  Single  variable  curves  may  be  generated  and 
studied,  e.  g.  capacity  as  a  function  of  temperature  at  a  constant  vibrational 
force  anywhere  from  10  to  kO  g. 

For  the  present  purpose,  however,  it  is  now  possible  to  make  the  required 
prediction  of  the  battery  design  reliability  at  the  temperature  requirement 
of  165  °F  and  at  the  vibration  requirement  of  10  g.  The  reliability  standard, 
it  is  estimated,  will  be  met  anywhere  to  the  left  of  the  reliability  boundary. 
Therefore  the  design  has  the  desired  reliability  at  the  l65°F,  10  g  point. 

If  desired,  it  is  possible  to  determine  how  many  standard  error  units  there 
are  along  the  Y  axis  through  (-1,  -l)  between  the  mean  response  surface  and 
the  10  second  plane  and  from  this  give  an  estimate  of  the  actual  reliability, 
99*999  etc.  However,  this  is  of  no  practical  significance,  the  main  question 
being,  has  the  reliability  standard  been  met  or  not. 

From  the  reliability  boundary  plot,  estimates  can  be  made  of  the  battery 
design  capability  at  environmental  stress  levels  above  the  requirements. 

These  estimates  are  qualified,  of  course,  by  the  reliability  standards.  Thus, 
expressed  with  95#  confidence,  it  may  be  estimated  that  no  more  than  one 
battery  out  of  a  thousand  will  fail  to  give  at  least  10  seconds  of  service 
under  the  following  environmental  conditions: 

a.  At  a  temperature  of  up  to  l88*F  while  undergoing  a  10  g 
vibration 

b.  In  excess  of  a  1(0  g  vibrational  force  at  a  temperature  of  165 *F 

c.  At  a  temperature  of  up  to  1®L°F  while  undergoing  a  20  g 
vibration,  etc. 
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Derivatio^of^xgerlEegte^^est^ata^ 

In  the  detailed  example  just  completed,  tests  vere  conducted  throughout 
the  X^,  space,  the  data  was  analyzed  and  various  estimates  were  made 
concerning  the  capabilities  of  the  population  from  which  the  sample  was  drown. 
These  estimates  in  final  form  are  represented  by  the  reliability  boundary 
contour.  A  logical  question  at  this  point  would  be,  how  close  is  this 
contour  to  the  true  picture?  For  example,  how  close  is  it  to  the  average 
contour  for  an  infinite  number  of  samples  of  the  same  size  drown  from  the 
same  population  (in  this  case,  future  batteries  produced  exactly  the  same 
in  every  respect  as  the  let  from  which  the  test  sample  was  drown)? 

There  are  two  ways  to  answer  such  a  question.  One  would  be  to  produce 
enough  batteries  and  run  enough  tests  with  samples  of  size  K  *  30  (l4  at 
(-1,  -l)  and  2  at  each  of  the  other  8  points)  to  begin  to  understand  the 
variability  of  the  individual  sample  reliability  boundaries .  In  this  way, 
it  could  eventually  be  said  whether  or  not  the  sample  gave  a  "good"  estimate 
of  the  true  picture.  This  approach  is  not  practical,  of  course,  with 
expensive  items  and  testing.  The  second  way  to  answer  the  question  is  to 
know  in  advance  the  nature  of  the  true  picture,  by  conducting  experiments 
with  data  which  is  derived  from  populations  of  known  statistical  parameters 
such  as  mean,  standard  deviation  and,  in  this  case,  coefficients  for  the 
two  variable  second  order  regression,  and  the  standard  error.  This  will  make 
it  possible  to  conduct  many  relatively  inexpensive  experiments  which  can  be 
made  to  closely  simulate  tests  with  actual  one-shot  items.  The  use  of  this 
technique  has  made  it  possible  to  thoroughly  explore  the  RSD  method  in  the 
Vital  area  of  the  effect  of  sample  size.  Completely  unexpected  phenomena 
were  observed  in  this  manner,  making  it  possible  to  derive  methods  for 
handling  them. 

In  establishing  the  simulated  test  data,  a  hypothetical  battery  design 
was  visualized  with  a  mean  capacity  of  l6.04  seconds  and  a  standard  deviation 
of  1.0  second  at  the  (-1,  -l)  requirement  point.  The  capacity  Y  of  the 
battery  was  established  as  a  function  of  X^  (Temperature)  and  X2  (Vibration 
g  Force)  in  accordance  with  the  following  equation: 

Y  -  14.61  -  3.24Xx  -  I.69X2  -  2.13X-12  -  O^Xg2  -  0.43X3X2 

A  variance,  standard  error  squared,  of  1.0  was  selected  for  the  X. ,  Xg  test 
region.  For  N  =  30  with  24  degrees  of  freedom  for  establishing  the  pooled 
standard  error,  the  K  factor  is  4.17.  Therefore,  KS  is  4.17,  and  a 
reliability  boundary  equation  may  be  set  up: 

3.2hXx  +  1.69X2  +  2.13X32+  0.94X3^  +  0.43X3X2  =  l4.6l  -  10.00  -  4.17  -  0.44 

The  contour  generated  from  this  equation  is  shown  on  Figure  9  as  (B),  the 
average  reliability  boundary  for  an  infinite  number  of  samples  of  size  N  *  30, 
d.f .  -  24.  It  is  by  no  means  the  reliability  boundary  for  an  infinite  sample 
size.  This  limiting  contour  is  generated  from  an  equation  similar  to  the 
above  except  that  KB  =  3.0 9,  where  the  K  factor  of  3*09  is  the  limiting 
value  of  K  for  an  infinite  sample  size  with  an  infinite  number  of  degrees  of 
freedom.  This  infinite  sample  size  contour,  or  universe  reliability  boundary, 
is  discussed  and  presented  later  in  Figure  17.  Curve  (a)  of  Figure  9  is  the 
mean  failure  contour  for  any  number  of  samples  of  any  size*  This  contour  is 
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generated  by  setting  Y  «  10  in  the  basic  equation,  giving: 

3.243^  +  1.69X3  +  2.13x3s  +  O^Xg2  +  0.43X3X3  -  l4.6l  -  10  -  4.6l 

It  is  seen  on  Figure  9  that  the  sample  mean  failure  contour  and  reliability 
boundary  obtained  in  Experiment  1  are  relatively  good  approximations  of  the 
true  values  representing  the  known  population.  Just  how  good  they  are  is  a 
question  that  could  best  be  answered  following  the  experimental  exploration 
of  the  effect  of  sample  size  in  determining  the  response  surfaces. 

Curve  (A)  requires  some  Comments.  It  passes  through  the  X^  axis  just 
about  at  the  limit  of  the  test  area.  It  passes  through  the  X3  axis  at  some 
point  well  beyond  the  test  area.  The  coefficients  of  the  equation  which 
generates  (a)  were  chosen  deliberately  So  that  the  curve  would  have  this 
general  shape.  This  represents  practical  considerations  in  establishing 
the  test  area,  as  previously  explained,  with  limits  either  at  the  estimated 
mean  failure  point  or  at  the  TEC  (test  equipment  capability)  • 


The  method  of  obtaining  actual  simulated  test  data  is  as  follows: 


The  base  equation,  Y  -  l4.6l  -  3.24X1  -  1.69X3  -  2.13X12  -  0.94Xg2 
0.43X3X3,  will  yield  the  following  predicted  values  of  Y  when  coded  3^  and 
Xj>  values  of  -1,  0,  or  1  are  substituted: 


a  table  of  random  normal  numbers5 with  a  universe  mean  of  0.00  and  a 
universe  standard  deviation  of  1.00  was  entered  at  random  in  order  to  derive 
the  test  data.  In  Example  1,  for  example,  colimn  Ho.  33  was  entered.  The 
first  17  numbers  are: 


,0.30 

-0.57 

1.62 

,0.43 

•0.41 

-1.53 

0.39 

-1.71 

•1.63 

2.75 

0.34 

1.48 

0.57 

,1.64 

0.70 

0.00 

I.69 

These  numbers  added  to  l6.04  give  the  following  Y  responses  at  (-1,  -l): 

15.74  15.47  17.66  15.61 

15.63  14.51  16.43  14.33 

14.41  18.79  16.38  17.52 

l6.6l  14.40  16.74  16.04 

17.73 
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These  numbers  have  a  mean  value  of  16.12  and  a  standard  deviation  of  1.313 
•which  are  estimates  of  the  true  population  values  of  l6.04  and  1.000. 

With  an  actual  large  lot  of  batteries  with  a  mean  of  16.04  and  and  a 
Standard  deviation  of  1*000,  it  is  possible  that  a  sample  of  N  =  17  could 
have  been  drawn  with  a  Y  of  16*12  and  an  S  of  1.313.  In  other  words,  the 
artificial  populations  can  be  as  useful  as  actual  populations  in  studying 
the  effects  of  performance  variability,  actually  much  more  so  since  all 
types  of  simulated  test  data  in  almost  unlimited  quantities  can  be  rapidly 
and  economically  generated. 

It  was  decided  to  explore  the  effects  of  Sample  size  in  using  the  RSD 
method  by  constructing  samples  with  three  representative  sanple  sizes: 

N  =  15,  N  =  30  and  N  =  60,  with  (-1,  -l)  point  test  quantities  of  7,  l4 
and  28  respectively  and  with  1,  2  or  4  tests  at  each  of  the  other  eight 
points  of  the  3^  test  design.  Ten  examples  were  constructed  for  each  of 
the  three  sample  sizes  for  a  total  of  30  examples .  The  examples  were 
constructed  in  such  a  way  that  the  results  of  the  analysis  could  be 
followed  in  stages  for  each  of  the  10  examples,  i.e,  from  N  e  15  to  N  =  30 
and  N  =  60  with  each  increase  based  on  additions  to  the  previous  data. 

A  table  of  random  numbers5 was  entered  to  give  10  consecutive  numbers 
of  6l,  46,  10,  24,  85,  ho,  38,  28,  58  and  17.  These  numbers  were  used  to 
determine  the  column  numbers  to  be  used  from  the  previously  referenced 
table  of  random  normal  numbers.  Thus  the  three  No.  1  samples,  1A  (N  =15), 
IB  (N  =  30),  and  1C  (N  =  60)  were  constructed  with  the  50  numbers  of  column 
6l  and  the  first  10  of  column  62.  The  No.  2  samples  were  derived  from 
columns  46  and  47,  etc. 

The  breakdown  of  the  60  random  normal  numbers,  e.g.,  from  columns  6l 
and  62,  was  as  follows,  with  the  random  numbers  added  to  the  predicted  value 
for  each  of  the  nine  test  design  points: 


Test  Design  Point 

Sample  Numbers 

%,  Xg 

1A  to  10A 

U2E2Es^! 

(-1,  -!) 

1-7 

mgm 

mam 

(  0,  -1) 

29 

h  rl) 

.  33 

■lltS’J 

“  (-1,  oy 

37 

(  Q,  0) 

4l 

mmm 

(  1,  Q) 

45 

IHEaCaH 

Ui,  1) 

49 

mmjsfwgm 

(  0,  1) 

53 

KiK 

HSUSE 

(  1,  1) 

57 

BB' 

As  an  example  of  the  above,  the  detailed  Example  1  was  based  on  the  7A 
and  7B  samples.  The  No.  7  samples  were  derived  from  columns  38  and  39  of 
the  table  of  random  normal  numbers.  Taking  the  (0,  0)  point  as  an  example, 
the  predicted  value  from  the  basic  equation  is  l4.6l.  Numbers  4l  and  42  in 
column  38  are  -1.75  and  +1.23.  Adding  these  to  l4,6l  gives  the  simulated 
test  data  values  of  12.36  and  15.84,  the  (0,  0)  values  used  in  the  example. 
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Some  of  the  factors  which  influence  sample  size  were  developed  in  the 
explanation  of  Example  1.  It  was  seen  that  testing  of  only  seven  units  at 
the  temperature  and  vibration  requirement  point  ( -1,  *1)  gave  an  erroneous 
indication  of  the  reliability  of  the  design  and  that  it  Was  necessary  to  test 
an  additional  16  units  before  the  desired  reliability  was  indicated. 
Similarly,  testing  only  one  unit  at  each  of  the  other  eight  points  of  the 
test  design  gave  an  erroneous  indication  of  the  shape  of  the  response 
surface,  implying  that  the  quadratic  components  were  not  significant. 

It  is  desirable  to  keep  the  required  sample  size  as  small  as  possible, 
particularly  when  dealing  with  expensive  units  and  tests .  In  order  to  gain 
insight  into  the  minimum  sample  size  which  will  give  effective  results, 
extensive  empirical  data  is  required,  the  purpose  of  this  section.  Thirty 
samples  were  established  in  accordance  with  the  procedures  given  in  the 
last  section.  The  three  sample  sizes,  15,  30  and  60,  were  arbitrarily 
chosen  as  practical  round  numbers .  The  established  test  data  for  the 
samples  is  given  in  Appendix  IIA.  The  remainder  of  Appendix  II  gives 
calculations ,  similar  to  those  in  Appendix  I  for  Example  1,  for  the  complete 
analysis  of  the  test  data  for  the  30  samples,  ending  with  the  reliability 
boundary  equations  in  Appendices  IH)6,  IIE6  and  IIP 5  from  which  the  contours 
were  plotted  as  shown  in  Figures  14,  15  and  16. 

The  inadequacy  of  a  subsample  of  7  units  at  (-1,  -l)  is  clearly  shown 
in  Appendix  IIB  where  5  of  the  10  "A"  samples  show  inadequate  reliability, 
i.e.  Y  -  KS  values  below  the  10  second  requirement.  This  occurred  with 
only  one  of  the  "B"  samples  (7B)  with  lU  tests  at  (-1,  *1).  As  shown  in  the 
previous  Example  1,  based  on  the  7A  and  7B  samples,  three  more  tests  were 
required  to  establish  the  reliability.  All  of  the  "C"  samples  with  28  tests 
at  (-1,  -l)  were  adequate  in  this  respect. 

The  inadequacy  of  the  sub  sample  of  7  at  (-1,  -l)  is  further  illustrated 
by  the  tests  of  the  assumption  of  normal  distributions  as  shown  on  Figures 
10,  11,  12  and  13  based  on  calculations  given  in  Appendix  UC.  Figure  10 
shows  in  general  a  wide  variation  of  the  cumulative  distribution  points 
around  the  best,  or  least  squares,  straight  lines,  indicating  little 
confidence  in  making  the  assumption  of  normality  even  when  it  is  known,  as 
in  this  case,  that  the  points  were  derived  from  a  known  normal  distribution* 
Figure  11,  showing  the  plots  for  the  "B"  subsamples  of  14,  indicates  that 
this  number  is  about  sufficient  for  making  the  required  assumption  of 
normality.  Most  of  the  solid  points  are  close  to  the  required  straight 
lines.  The  up,  down  and  up  plot  of  sample  5B  is  characteristic  of  a  double 
population  distribution  on  either  side  of  the  mean  value.  This  actually 
occurs  in  the  column  of  random  normal  numbers  from  which  the  sample  was 
derived.  In  practice,  however,  there  is  no  known  way  in  which  a  battery 
performance  parameter  can  be  distributed  in  this  manner  other  than  through 
chance  sample  variability.  Therefore  the  plot  of  5B  could  be  considered  to 
represent  a  straight  line,  thus  satisfying  the  assumption  of  normality. 

The  tendency  most  to  be  looked  for  in  analyzing  these  plots  is  a  pronounced 
curvature  in  one  direction  or  the  other,  suggesting  the  probability  of  a 
skewed  distribution.  A  slight  tendency  of  this  sort  is  noted  in  the  plot 
of  1QB,  but  not  enough  to  counter  the  assumption  of  normality.  The  plots 
for  the  "C"  samples  in  Figures  12  and  13  all  tend  to  satisfy  the  assunptio- 
of  normality  of  distribution. 


Figure  12 

28 


Figure  13 
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All  of  the  cumulative  distribution  plots  have  dotted  lines  on  them,  all 
parallel  to  each  other.  These  represent  the  straight  line  plots  for  an 
infinite  sample  size  with  a  Standard  deviation  of  One  and  are  shown  to 
indicate  how  veil  the  individual  samples  Correspond  to  the  true  mean  and 
standard  deviation.  It  is  seen  that  the  larger  the  sample  the  closer  the 
correspondence  in  general.  The  three  Vertical  lines  in  each  figure  represent 
the  mean,  at  P  »  50,  and  the  minus  and  plus  one  standard  deviation  values  at 
P  «  16  and  P  -  84. 

Based  on  the  foregoing,  it  may  he  stated  that  when  dealing  with  samples 
drawn  from  a  normally  distributed  population  with  relatively  low  variance, 
a  sample  of  approximately  14  at  (-1,  -l)  will  in  general  Satisfy  the 
assumption  of  normality  with  a  cumulative  distribution  plot  which  approximates 
a  straight  line.  Slight  departures  from  normality  cannot  reasonably  be 
detected,  but  the  probability  is  that  relatively  large  degrees  of  skewness 
can  be.  Much  more  confidence  in  assuming  a  normal  distribution  can  of  course 
be  had  with  a  larger  sample  size  if  this  is  economically  feasible. 

In  conducting  act  actual  experiment,  once  it  has  been  determined  that  it 
is  not  unreasonable  to  assume  normality  of  distribution,  at  least  in  the 
region  about  (-1,  -l),  the  remainder  of  the  test  region  will  be  explored  to 
estimate  the  nature  of  the  response  Surface.  The  minimum  subsample  size  for 
this  is  eight,  one  at  each  of  the  remaining  test  design  points.  Even  for  a 
sample  drawn  from  a  population  with  relatively  low  variance.  Example  1  showed 
that  this  minimum  subsample  size  is  inadequate.  This  IS  more  positively 
demonstrated  by  Figure  14,  the  "A"  sample  reliability  boundaries  drawn  from 
the  calculations  given  in  Appendix  HD.  Nine  of  the  10  boundaries  fall 
above  (-1,  -l)  but  the  wide  variation  of  them  around  curve  B  Indicates  a  high 
probability  of  one  of  these  small  samples  indicating  that  the  battery  design 
is  not  reliable  when  it  actually  is.  Sample  7A  is  in  this  category.  The 
mean  failure  contour  for  the  sample,  as  shown  by  the  dots,  gives  a  good 
estimate  of  the  true  mean  failure  contour  A.  In  fact  each  of  the  10  samples 
do,  as  shown  by  the  cluster  of  points  around  curve  A  at  Xg  =  0.  The  wide 
variation  in  the  reliability  boundaries  is  due  to  the  variation  which  the 
samples  give  in  the  estimate  of  the  standard  error.  There  is  also  a  wide 
variation  in  the  estimates  of  the  coefficients  of  the  response  surface 
equations  (Appendix  HDl) .  This  accounts  for  the  considerable  variation 
in  the  shapes  of  the  contours.  Another  indication  of  the  Inadequacy  of  the 
sample  size  is  given  in  the  Analysis  of  Variance  in  Appendix  HD2.  Only 
four  of  the  samples  show  the  quadratic  terms  to  be  significant.  For  the 
others  linear  response  surfaces,  or  planes,  are  erroneously  indicated  with 
straight  line  contours. 

Replication,  repeating  the  experiment,  is  found  to  eliminate  all  of  these 
problems  as  shown  with  the  reliability  boundaries  for  the  "B"  samples  as 
given  in  Figure  15,  based  on  the  calculations  in  Appendix  HE.  It  is  seen 
that  each  of  the  10  samples  now  gives  a  fairly  good  estimate  of  the  true 
contour,  represented  by  curve  B.  The  Analysis  of  Variance  given  in  Appendix 
HE2  shows  that  each  sample  correctly  Indicates  that  the  quadratic  components 
of  the  mathematical  model  are  significant.  Having  two  results  at  each  test 
design  point  also  affords  a  measure  of  the  fit  of  the  mathematical  model  to 
the  data.  In  every  case  the  lack  of  fit  term  is  found  to  he  Insignificant. 

The  assumption  of  uniformity  of  variance  throughout  the  entire  test  region 
Is  satisfactorily  met  for  each  sample  as  shown  in  Appendix  XXE5.  Again, 
uniform  variance  is  inferred  by  the  Inability  to  prove  the  variance  at  (-1,-1) 
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A*  Universe  Mean  Failure  Contour 
B*Average  R.B.,  Infinite  Number  of  Samples 


Figure  14 
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IB  to  10B  Sample  Reliability  Boundaries, 2 Tests  per  Design  Point 


to  be  significantly  different  from  the  variance  throughout  the  rest  of  the 
test  region.  This  mates  it  possible  to  compute  a  uniform  pooled  standard 
error  based  on  all  of  the  test  data  and  from  this  to  calculate  the  equations 
for  the  reliability  boundaries. 

The  reliability  boundaries  for  the  "C"  samples,  with  four  tests  at  each 
of  the  eight  test  points,  as  given  in  Figure  16,  naturally  show  a  still 
Closer  estimate  of  the  true  curve  B.  Whether  they  give  that  much  better  an 
estimate  than  the  "B"  samples  do  to  justify  the  added  costs  is  a  debatable 
point.  In  the  opinion  of  the  author,  they  do  not.  However,  for  a  battery 
design  with  much  greater  variability  than  the  hypothetical  design,  or  for 
samples  which  give  reliability  boundaries  on  or  slightly  below  (-1,  -l),  it 
would  undoubtedly  be  necessary  to  replicate  beyond  two  tests  at  each  point. 

An  analysis  of  variance  is  not  included  in  Appendix  HF  for  the  "C" 
samples  since  the  analysis  for  the  "B"  samples  had  shewn  these  to  be 
Satisfactory  in  regard  to  significance  of  the  quadratic  terms  and  proper 
fit  of  the  mathematical  model. 


The  effects  of  replication  are  shown  quite  clearly  on  Figure  17.  The 
"A"  sample  reliability  boundary  ?.l  gives  &  poor  estimate  of  the  true  curve 
1A  (curve  B  of  Figure  lU).  The  "B"  sample  boundary  gives  a  much  better 
estimate  of  the  true  curve  2A  and  this  is  even  more  pronounced  with  the  "C" 
sample  curves  and  kA.  Additional  curves  are  Shown  for  a  tremendous  sample 
consisting  of  40  tests  at  each  of  the  eight  design  points.  Contour  40 
corresponds  very  closely,  as  one  would  expect,  with  !+OA  representing  the 
average  reliability  boundary  of  an  infinite  number  of  samples  with  Uo  at 
each  point.  Curve  ko  M.F.C.,  the  mean  failure  contour.  Is  found  to 
correspond  perfectly  with  the  total  population  or  universe  mean  failure 
contour*  The  curves  1A,  2A,  4a,  and  kOA  are  seen  to  approach  a  limiting 
contour  which  is  called  the  U.R.B.,  or  universe  reliability  boundary.  This 
curve  is  derived  from  an  equation  based  on  a  K  factor  of  3*09  for  a  sample 
size  of  Infinity.  It  may  be  interpreted  by  saying  that  in  testing  an 
infinite  sample  size,  one  out  of  a  thousand  units  would  fall  below  the  U.B.B. 
Fifty  per  cent  would,  of  course,  fall  below  the  universe  mean  failure  contour. 
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Tunpyrotutf,  °F 


Effect  of  Replication  on  R.B. 


CONCLUSIONS 

The  theory  of  the  PSD  method  has  "been  covered,  the  mathematics  involved 
have  been  explained  in  detail,  and  an  analysis  of  the  effects  of  Sample  Size 
has  been  performed  in  order  to  provide  a  guide  in  determining  sample  Sizes 
for  an  actual  reliability  program  for  batteries  or  other  one-shot  items.  It 
has  been  demonstrated  that  the  method  will  provide  a  great  deal  of  valuable 
reliability  prediction  data  from  the  testing  of  a  relatively  small  sample. 

As  posed  at  the  end  of  the  first  section  of  the  discussion,  a  basic 
question  is,  what  assurance  is  there  that  the  final  reliability  predictions 
are  correct  and  usable?  This  is  particularly  relevant  since  the  sample  test 
data  was  not  only  used  for  making  the  predictions,  but  for  testing  the 
assumptions  of  normality  of  distribution  and  uniformity  of  variance  upon 
which  basis  the  predictions  were  made.  The  sample  size  determining  section 
showed  that  the  sample  size  has  great  bearing  on  the  assurance  with  which 
the  final  conclusions  may  be  accepted.  The  smallest  sample  size  of  15  showed 
a  relatively  high  probability  of  erroneously  concluding  that  the  battery 
design  was  unreliable .  The  wide  scatter  of  the  sample  reliability  boundaries 
also  showed  that  little  faith  could  be  placed  in  reliability  estimates  out 
along  the  X  axes.  The  predictions  for  N  =  60  showed  such  excellent  conform¬ 
ance  with  the  true  picture  that  the  reliability  predictions  may  be  made  with 
great  assurance  in  their  accuracy.  However,  at  least  for  most  missile 
battery  programs,  a  Sample  size  of  60  would  be  Considered  excessive  from  an 
economic  standpoint.  Forty-five  would  probably  be  as  well.  Thirty  would  be 
considered  a  reasonable  size  in  many  caseB. 

The  results  for  the  "B"  samples  of  size  30  showed  in  general  a  good 
approximation  of  the  characteristics  of  the  battery  population.  The  test  of 
the  suitablility  of  assuming  normal  distribution  appeared  satisfactory,  as 
well  as  the  test  for  significance  of  quadratic  components,  the  test  of  the 
assumption  of  unifora  variance  and  the  actual  reliability  boundaries 
themselves.  These  factors  indicate  that  the  sample  size  of  30  should  be 
adequate  in  most  cases  and  that  the  final  reliability  predictions  can  be 
with  considerable  assurance.  In  conducting  an  actual  program,  as  previously 
pointed  out,  there  may  be  instances  where  more  than  14  tests  will  be  required 
at  the  Xj,  Xg  requirement  point  and  where  it  may  be  necessary  to  replicate 
beyond  two  tests  per  point  throughout  the  rest  of  the  test  design.  In 
establishing  the  total  lot  size  for  the  program,  a  reserve  should  be  provided 
for  these  contingencies . 
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AHESDIX  I 


Exgeriment^_ 


B.  Cumulative  distribution  points  of  responses  at  X 


Select  suitable  intervals  for  X  such  that  there  are  from  10  to  20 
intervals.  Then  set  up  a  table  for  the  midpoints  of  the  selected  T 
intervals,  tally  the  T  responses  and  set  up  three  additional  columns  for 
frequency,  cumulative  frequency,  and  cumulative  frequency  expressed  as 
percentages  of  2  x  N  *  3 h.  The  cumulative  frequency  column  is  derived  from 
the  frequency  column  by  a  counterclockwise  rotating  addition  system;  e.g.. 


0+2 
+ 

2«U,  etc. 

The  cumulative  frequency  values  thus  obtained  are  divided  by  the  double  total 
to  yield  the  cumulative  frequency  percentages  (expressed  to  tenths  of  a  per 
cent  below  1 0%  end  above  9050* 
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XX-.T  -22,00 
*1-—  "  — - 

E3UT  -10.03 

Ij.  m  b  as - r— i rr 

2  6  6 

XXjXjT  -U.50 


12 


bll“ 


b22  " 


'  -1.13 

U  U 

XxA-2/3Tt  75.78-77.29 


2 XA  T-2/32T  76, 89-77.29 

2 


-1,51 

2 

-0.1(0 

2 


-0.76 


-0.20 


b0-^|  -  2/3bn-2/3b22  -  12.881  +  0.507  ♦  0.133  -  13. 


Therefore  the  fitted  equation  is: 


T  -  13.52-3.67X1-1.67X2-0,76X12-0.2aX22-1.13X1X2 


A  check  is  made  for  arithmetical  errors.  The  sum  of  the  coefficients  is  6.09. 
This  is  checked  against  the  addition  of  the  X  response  values,  ordered  from 
X^  to  ,  in  accordance  with  the  following  formula: 

5\^(VVV-V8(Vt8)-V29I9 

36 


5*l6.12-h(l6.07+l6.99+12.86)-11.27+8(9.19+11.22)-2S.7&*29!x6.1(3 


80.60-183. 68-11.27+163. 28-15. 78+186.U7  _  219.62  _  ,  _ 

*  - - - _ .  1 - - -  ■  • -  *  — -  ■  0.10 
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Allowing  for  a  difference  due  to  rounding  off  errors,  this  cheek  indicates 
that  no  arithmetical  errors  hams  been  ^ade  In  aatabliwhing  the  coefficients' 


0.9U21 


2.0 9h2 

F,  ,  -  - ; - -  2.22 

3,3  o.9h21 


In  order  for  the  quadratic  components  to  be  significant  at  a  95%  significance 
level,  the  F_  -  ratio  most  be  at  least  9.28.  The  ratio  of  2,22  shows  the 
quadratic  components  to  be  not  significant.  Another  way  of  stating 

this  is  that  the  quadratic  components  hare  not  been  proved  to  be  significantly 
different  from  zero. 


S.  Determination  of  coefficients  of  a  fitted  model  of  the  form? 


1  - v v  *v»*hA*v»* va 

hi 


ZX3I 

—  ■!■■■  at 

UU.29 

■■  ■■■  ■, ■■■■-.  m 

-  3.69 

12 

IXgT 

12 

20*31 

_  i -Ao 

12 

12 

IXjXgl 

7.22 

-  -  O.J 

8 

8 

b  _  ■ 

n 


22 


IX,  ^-2/311  lU7.li9-35U.05  -6.56 

U  ’  ~ 

-15U.Q5  2.36 


«■  -  i*6U 


-  0.59 


II 

b  -  —  -  2/3b  -2/3b  -  12.838  +  1.093  +  0.393  -  1U.32 

0  x8  11  22 


Therefore  the  fitted  equation  is: 

T  -  1U.32  -  3,691!  -  1.6912  -  l.^2  -  -  0.901^2 


Cheek:  ^Coefficients  *5.81 
Coefficients  » 


53^  -^VW  -Y3  +8(Y6-Hq)  -X?  +29Yy 
72  (for  N  -  18) 


5(32.2U)  -U(9U.10)  -21.1U  +8(li0.18)  -30.87  +29(12 .55) 

72 

161.20  -376.UO  -21.1U  +321.1iU  -30.87  +363.95 
72 


5.81 


U3 


U16.18 

72 


-1 


The  critical  F,  ratio  at  a  95%  level  of  significance  is  3*1*9.  This  is 
exceeded  by  the  6.50  ratio.  Therefore,  the  quadratic  components  are 

significant.  The  lack  of  fit  ratio,  1.U5,  is  less  than  Pi  g  »  3.8U,  there¬ 
fore  the  lack  of  fit  term  is  not  significantly  different  u*  fTon  experi¬ 
mental  error  and,  therefore,  is  non-significant. 


-  1.00893 


5  d.f .  are  expended  for  estimating 
the  coefficients;  therefore,  the 
denominator  is  N-f>  *  16-5  ■  11  . 
The  sixth  coefficient  is  based  on 
a  d.f.  iron  the  tests  at  (-1,  -1). 
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H. ...  Standard  error  from  pooling  of  estimates  of  yagianee 


ll(d.f.)  x  1.00893  *  l6(d.f.)  x 
27  d.f. 

11.09823  27.59296  _  38.69119 

27  27 


Sp  (pooled  standard  error)  • 


*  1.197 


U6 


L.72U56 


-  1.13301 


I.  Determination  of  Mean  Failure  Contour  Points 


1.  Mean  failure  contour  equation: 


3.70^  +  1.70X2  +  1.65XJ2  +  O.toXg2  +  0.92Xj,Xg  *  U.33 

2.  Value  of  X^  at  Xg  -  0  ;  3.703^  +  1.651^  •  U.33 

estimate  X^  at  0.80  ;  2.96  +  1.06  •  U.02 

estimate  X^  at  0.85  j  3.15  +  1.19  *  U.3U 

estimate  X^  at  0.8U  $  3»ll  +  1.16  •  U.27  ,  •*.  root  is  0.85  *  ^ 

3.  Value  of  Xg  at  X^  •  0  ;  1.70X2  +  0.60X22  «  U.33 

estimate  X2  at  2.5  (beyond  test  region,  maximum  X2  »  1.0) 

U.  Value  of  X1  at  X2  »  1  j  3.700^  +  1.70  +  1.6£X12  ♦  0.60  ♦  0.923^ 

■  U.33 

U.62X1  +  1.65X12  -  U.33  -  1.70  -  0.60 


estimate  X^  at  0.U0  j  1.85  +  0.26  »  2.11 

estimate  X1  at  O.38  ;  1.76  +  0.2U  ■  2.00 

estimate  X^  at  0.39  j  1.80  +  0.25  -  2.Q5  $  .*.  root  is  0.39  -  X^ 

5.  Value  of  X1  at  Xg  «  -1  j  2.780^  +  X.65^2  -  U.33  +  1.70  -  0.60 
■  5*1*3 

estimate  X^  at  1.50  j  U.17  +  3.71  -  7.88 

estimate  at  1.20  j  3.3U  +  2.38  -  5.72 

estimate  X^  at  1.10  j  3.06  +  2.00  -  5.06 

estimate  ^  at  1.15  ;  3*20  +  2.18  -  5.38 

estimate  X^  at  1.16  ;  3.22  +  2.22  -  5.1»U  ,  .*.  root  is  1.16  - 

MSM L  Extraction  of  the  roots  by  the  above  trial-and-error  method  is  very 
fast  with  a  desk  calculator,  much  faster  than  using  the  standard 
formula  JSTrT 


DeterBdnationofreliabillt^boupda^fY^^^eontomjgoints^ 


X.  Reliability  boundary  equation; 

3.70^  +  1.70X2  +  +  0.60Xg2  +  0.92XjXg  ■  -  0.57 

2.  Value  of  2^  at  Xg  •  0  ;  3.7CX1  +  1.65X12  -  -  0.57 

estimate  X^  at  -0.20  ;  -  0.7U  +  0.07  *  -  0.67 

estimate  X^  at  *0.18  ;  -  0.67  +  0*05  *  -  0.62 

estimate  X^  at  -0.16  ;  -  0*59  +  0.0U  *  -  0.55 

estimate  X^  at  -0.1?  ;  -  0.63  +  0.05  8  -  0.58  ,  .*.  root  is 

-0.17  »  Xj^ 

3.  Value  of  Xg  at  Xj^  -  0  ;  1.70X2  +  0.60Xg2  -  -  0.57 
estimate  Xg  at  -0.U0  j  -  0.68  +  0.10  -  -  0.58 
estimate  Xg  at  -0*39  5  -0.66  +  0.09  8  -  0.57  »  root  is 


U. 


5. 


6. 


-  0.39  8  Xg 


2 

Value  of  X^  at  Xg  -  -1  ;  2.781^^  +  1.6$!^  -  -  0.57  +  1.70 

~ ~ -  0.60  •  0.53 
estimate  X^  at  0.20  ;  0.56  +  0.07  -  O.63 

estimate  X^  at  0.18  ;  0.50  +  0.05  -  0.55 

estimate  X^  at  0.17  ;  0.U7  +  0.05  -  0.52  ,  .*.  root  is  0.17  8  X^ 

Value  of  \  at  Xg  -  0.5  5  U.l6xi  +  1.65XL2  «*  -  0.57  -  0.85^_ 

-  2.08  +  0.U1  -  -  1.67 
-1.66  +  0.26  -  -  1.50 

-  1.87  ♦  0.33  8  -  1.5b 

-  1.96  +  0.36  8  -  1.60  ,  root  is 

-  O.U6  -  X^^ 


estimate  X^  at  -  0.50  j 
estimate  X^  at  -  0.U0  j 
estimate  X^  at  -  0.b5  J 
estimate  X_  at  -  0.b7  j 


Value  of  X^  at  Xg  «  1.0  ;  U.62X^  +  1.65Xg  *  -  0.57 

estimate  X^  at  -1.0  j  -  U.62  +  1.65  8  -  2.97 
estimate  X^  at  -0.96  ;  -  b.bb  +  1.52  *  -  2.92 


estimate  X^  at  -0.9b  j  -U.3U  +  l.b6  *  -2.88 


-  1.70 
8  -  2.87 


estimate  X  at  -0.93  f  -  b.30  ♦  1.1(3 
1  b8 


2.87  t  ••  root  is 
-  0.93  8  Xx 


APPEHUEX  II 

A.  Sample  Test  Data 
1.  "A"  Samples 


S  AM  P  L  E 


Test 

Design 

Point 

1A 

2A 

E9 

mm 

ns 

ns 

ns 

8a 

mm 

ess 

(-1,-1) 

16.09 

15.72 
16.98 
16. 5U 
18.35 
17.86 
15.62 

15.21 

16.40 

16.54 

14.55 
14.99 
14.78 
17.06 

1 

1 

16.36 

14.95 

16.94 

17.07 

14.74 

15.38 

17.39 

14.48 

16.45 

17.28 

16.87 

16.99 

16.76 

14.76 

15.86 

15.87 
15.32 
17.45 

14.95 

15.85 

16.19 

15.74 

15.63 

14.41 

16.61 

15.47 

14.51 

18.79 

15.98 

15.50 

16.27 

15.96 

16.23 

14.83 

16.38 

15.55 

16.24 

16.22 

16.85 

16.56 
14.38 
14.37 

15.77 
16.06 
15. 04 
17.20 
15.54 
15.86 
14.79 

(0,-1) 

16.94 

16.44 

16.36 

15.18 

14.16 

16.53 

16.07 

l4.4o 

16.03 

ESS 

(1,-1) 

9.93 

n.30 

10.49 

II.05 

11.62 

9.21 

11.27 

11.72 

9.65 

m 

15.77 

15.80 

16.83 

14.93 

16.38 

16.32 

16.99 

15.88 

14.20 

14.62 

(0,0) 

13.40 

15.05 

14.31 

15.60 

14.50 

15.00 

12.86 

13.61 

13.51 

14.89 

(1,  0) 

7.71 

10.61 

8.04 

9.36 

8.13 

9.71 

9.19 

7.00 

8.03 

IE 

(-1,  1) 

13.77 

13.88 

12.37 

13.77 

13.40 

13.47 

15.78 

13.39 

12.75 

KTS«74 

V  Vnv 

(o,  1) 

11.92 

10.85 

12.60 

13.39 

11.18 

12.36 

11.22 

12.11 

12.43 

m 

(1, 1) 

6.61 

4.74 

5.61 

7.65 

6.66 

5.58 

6.43 

5.97 

4.14 

s 
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SAM  P  L  E 


2.  "B"  Samples 


Test 

D®  si  Eft 

Point 

IB 

2B 

3B 

4B 

5B 

6B 

7B 

8B 

2,.,,. 

1GB 

(-1,-1) 

16.09 

15.21 

15.08 

16.36 

14.48 

15.86 

15.74 

15.98 

15.55 

15.77 

15.72 

16.40 

16.57 

14.95 

16.45 

15.87 

15.63 

15.50 

16.24 

16.06 

16.98 

16.54 

14.17 

16.94 

17.28 

15.32 

14.41 

16.27 

16.22 

15.40 

16.5^ 

14.55 

15.77 

17.07 

16.87 

17.45 

16.61 

15.96 

16.85 

17.20 

18.35 

14.99 

16.00 

14.74 

16.99 

14.95 

15.47 

16.23 

16.56 

15.54 

17.86 

14.78 

16. 41 

15.38 

16.76 

15.85 

14.51 

14.83 

14.38 

15.86 

15.62 

17.06 

15.34 

17.39 

14.76 

16.19 

18.79 

16.38 

14.37 

14.79 

16.05 

17.36 

15.61 

15.61 

15.11 

17.32 

14.40 

16.86 

15.59 

14.94 

17.23 

15.68 

15.57 

16.05 

14.66 

18.50 

17.66 

17.17 

17.22 

16.35 

17*34 

15.18 

16.16 

13.53 

17.18 

16.56 

16.43 

15.28 

16.55 

16.25 

16.05 

15.84 

15.80 

15.89 

15.74 

15.59 

16.38 

16.68 

16.97 

15.51 

15*45 

14.88 

15.17 

15.91 

14.84 

16.84 

16.74 

14.60 

16.88 

15.21 

15.92 

15.73 

15.02 

15.43 

17.16 

14.40 

15.61 

16.04 

16.95 

16.95 

16.92 

13.90 

16. 46 

16.42 

15.77 

15.16 

14.33 

17.46 

15.84 

15.32 

(0,-1) 

16.94 

16.44 

16.36 

15.18 

l4.l6 

16.53 

16.07 

l4.4o 

16.03 

14.11 

14.74 

14.36 

14.34 

16.26 

13.95 

13.82 

16.55 

16.28 

14.56 

15.70 

(1,-1) 

9-93 

11.30 

10.49 

U.05 

11.62 

9.21 

11.27 

11.72 

9.65 

9.43 

10.20 

10.82 

10.94 

11.09 

9.72 

8.99 

9.87 

10.28 

12.29 

9.82 

(-1,0) 

15.77 

15.80 

16.83 

14.93 

16.38 

16.32 

16.99 

15.88 

14.20 

14.62 

14.83 

16.84 

16.79 

16.03 

15.29 

15.27 

15.79 

15.78 

14.48 

15.03 

(0,0) 

13.40 

15.05 

14.31 

15.60 

14.50 

15.00 

12.86 

13.61 

13.51 

14.89 

14.12 

16.54 

15.67 

13.59 

13.81 

14.39 

15.84 

12.69 

15.89 

14.28 

(1,0) 

7.71 

10.61 

8.04 

9.36 

8.13 

9.71 

9.19 

7.00 

8.03 

8.93 

9.14 

7.83 

10.34 

9.72 

7.57 

12.09 

8.72 

8.69 

11.62 

10.87 

(-1,1) 

13.77 

13.88 

12.37 

13.77 

13.40 

13.47 

15.78 

13.39 

12.75 

15.86 

13.54 

13.39 

14.82 

13.14 

13.90 

13.55 

15.09 

14.87 

15.60 

13.93 

(0,1) 

11.92 

10.85 

12.60 

13.39 

11.18 

12.36 

11.22 

12.11 

12.43 

11.06 

12.48 

11.24 

12.68 

13.71 

10.66 

13.45 

11.05 

9.62 

11.86 

13.22 

(1,1) 

6.61 

4.74 

5.61 

7.65 

6.66 

5.58 

6.43 

5.97 

4.14 

5.37 

6.74 

5,36 

6.28 

,  4.91 

..  »•* 

t  6.41 

6.12 

-  5’92 

5.27 

4.94 
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3*  "C"  Samples 


S  AM  PL  E 


I  Test 

1 

Point 

1C 

!■ 

sm 

Km 

■n 

6e 

m 

sm 

mm 

(-1,-1) 

16.09 

15.21 

15.08 

16.36 

14.48 

15.86 

15.74 

15.98 

15.55 

15.77 

15.72 

16.40 

16.57 

14.95 

16.45 

15.87 

15.63 

15.50 

16.24 

16.06 

16.98 

16.54 

14.17 

16.94 

17.28 

15.32 

l4.4l 

16.27 

16.22 

15. 4o 

16. 54 

14.55 

15.77 

17.07 

17.45 

16.61 

15.96 

16.85 

17.20 

18.35 

14.99 

16.00 

14.74 

16.99 

14.95 

15.47 

16.23 

16.56 

15.54 

17.86 

14.78 

16.41 

15.38 

16.76 

15.85 

14.51 

14.83 

14.38 

15.62 

17.06 

15.34 

17.39 

14.76 

16.19 

18.79 

16.38 

14.37 

14.79 

16.05 

17.36 

15.61 

15.61 

15.U 

17.32 

l4.4o 

16.86 

15.59 

14.94 

17.23 

15.68 

15.57 

16.05 

14.66 

18.50 

17.66 

17.17 

17.22 

16.35 

17.34 

15.18 

16.16 

13.53 

17.18 

16.56 

16.43 

15.28 

16.55 

16.25 

16.05 

15.84 

15.80 

15.89 

15.74 

15.59 

16.38 

16.68 

16.97 

15.51 

15.45 

14.88 

15.17 

15.91 

14.84 

16.84 

16.74 

14. 60 

16.88 

15.21 

15.92 

15.73 

15.02 

15.43 

17.16 

l4.4o 

15.61 

16. 04 

16.95 

16.95 

16.92 

13.90 

16.46 

16.42 

15.77 

15.16 

14.33 

17.46 

15.84 

15.32 

16.48 

14.26 

16.10 

16.43 

17.44 

15.76 

17.52 

15.50 

15.64 

15.34 

16.69 

15.26 

16.60 

16.57 

15.05 

14.85 

16.  o4 

17.28 

16.34 

16.65 

15.53 

15.00 

13.68 

16.38 

16.66 

15.62 

17.73 

15.37 

14.87 

15.30 

16.28 

15.53 

18.00 

14.69 

15.53 

14.43 

15.85 

16.48 

15.62 

16.09 

16.48 

15.98 

15.76 

17.16 

16.55 

15.06 

16.80 

16.87 

16.56 

17.81 

16.31 

16.21 

16.97 

16.75 

15.45 

16.53 

15.77 

17.34 

16.43 

17.09 

16.30 

16.23 

15.01 

14.87 

16.65 

14.92 

15.02 

14.97 

14.62 

14.85 

16.90 

14.87 

16.52 

16.30 

15.72 

16.69 

15.57 

15.57 

14.72 

16.80 

16. 14 

15.89 

18.75 

15.52 

15.67 

15.84 

16.53 

14.15 

14.58 

15.56 

16.56 

15.21 

15.98 

15.96 

16.45 

17.79 

16.95 

16.23 

17.51 

15.63 

15.73 

16.23 

15.74 

17.59 

16.61 

16.58 

16. 04 

16.50 

16.03 

14.70 

15.54 

16.91 

15.45 

16.63 

16.27 

14.49 

17.07 

16.73 

14.95 

15.13 

16.  o4 

15.68 

14.99 

16.38 

16.96 

15.52 

13.87 

16.74 

15.25 

15.19 

17.20 

15.62 

14.69 

16.26 

18.47 

15.82 

16.35 

15.99 

16.19 

17.52 

3.  "C"  Samples  (Coat'd) 


SAMPLE 


(0,-1) 


(1,-1) 


(0,0) 


(1,0) 


(-1,1) 


(0,1) 


(1,1) 


16.94  16.44  I6.36 
14.74  14.36  14.34 

16.94  15.78  15.91, 

15.99  15.34  15.78  15.16 


9.93  Hi  30  10.49 

10.20  10.82  10.94 
9.63  12.18  10.90 
10.99  10.51  10.75 


15.06 


7.71  10.61  8.04 
9.14  7.83  10.34 
9.31  8.77  8.32 
9.36  9.60  10.46 


13.77  13.88  12.37 
13.54  13.39  14.82 
15.02  15.14  12.19 
12.65  13.89  14.80 


11.62  9.21  11.27 
9.72  8.99  9.87 
9.54  10.87  9.10 
10.49  9.19  9.54 


14.94 


8.13  9.71 
7.57  12.09 
9.43  11.56 
9.85  8.03 


13.40  13.47 
13.90  13.55 
14.52  12.66 
12.07  13.24 


11.92  10.85  12.60  13.39  H.18 

12.48  11.24  12.68  13.71  10.66 

11.58  13.54  12.08  11.71  12.99 

12.08  12.07  10.60  H.39  12.16 


12.86  13.61 
15.84  12.69 
15.H  15.25 
14.22  1 16.46  14.45 


16.38 


12.29 

10.82 

IO.65 


14.20 

14.48 

17.00 

14.51 


8.03 

11.62 

8.07 

9.15 


8.93 

10.87 

10.49 

8.49 


13.39  12.75 
14.87  15.60 

13.94  13.20  1  13.20 
12.14  12.80  1 13.80 


11.22  12.11  12.43 
11.05  9.62  11.86 

11.33  10.90  12.45 
U.90  10.94  13.59 


C.  Testing  the  Assumption  of  Normal  I^tijbution_at_Polnt_X^^Ll.l 


The  necessary  calculations  are  given  for  Samples  1A,  IB  and  1G  for  7,  14, 
and  28  results  respectively  at  (-1,  -l).  Calculations  for  the  remaining 
samples  are  similar.  (See  Appendix  IB  for  specific  instructions). 


Sample  1A*  N 

*2 

Sample  ib 

,  N  =  14 

Samp; 

Le  1C,  N  m  28 

Y 

Freq 

Cum 

Free 

cum 

Freq 

of 

2N) 

Free  _ 

Cum 

Freq 

Cimn 

Freq 

of 

2N) 

Freq 

Cum 
„  Freq 

Cum 
Freq 
($  of 
2N) 

18.6 

0 

0 

0 

0 

_ 

0 

0 

18.2 

1 

1 

7.1 

n 

1 

3.6 

1 

1 

1.8 

17^8 

1 

3 

21 

WM 

3 

11 

1 

3 

5.4 

17.4 

0 

4 

29 

B3i 

6 

21 

3 

7 

13 

17.0 

1 

5 

36 

■Xi 

10 

36 

3 

13 

23 

1 6.6 

2 

8 

57 

■SB 

14 

50 

6 

22 

39 

16.2 

1 

IB 

79 

magi 

19 

68 

8 

36 

64 

15.8 

1 

92.9 

24 

86 

3 

47 

84 

15.4 

0 

mm 

• 

BflS 

27 

96.4 

3 

53 

94.6 

15.0 

mM 

0 

28 

- 

0 

56 

- 

The  cumulative  frequency  percentages  are  then  plotted  against  the  midpoint 


values  of  the  Y  intervals  on  arithmetical  probability  paper  as  shown  on 
Figures  10,  11,  12  and  13.  The  least  squares  straight  lines  are  then  drawn 
through  the  solid  points  representing  percentages  between  10  and  90-  In  the 
case  of  the  10  samples  1A  to  10A  on  Figure  10,  a  best  straight  line  can  be 
represented  only  for  sample  8A.  The  wide  divergence  of  the  points  for  the 
other  samples  require  calculation  of  the  equation  for  the  least  squares 
straight  line.  For  the  "B"  samples  this  had  to  be  done  for  samples  5B,  6b, 
7B,  9®  and  10B.  For  the  "C"  samples  it  had  to  be  done  only  for  5C  and  IOC, 
the  other  samples  permitting  the  estimation  of  the  best  straight  line 
directly.  A  technique  for  deriving  the  least  squares  straight  line  is 
illustrated  for  sample  5C. 

The  solid  points  for  the  sample  are  transferred,  from  Figure  12,  to 
regular  graph  paper. 


Arbitrary  scale  values  are  assigned  to  the  two  axes  called  X  and  Y  for 
convenience*  Values  of  X  and  Y  for  the  points  are  read  off  the  graph* 
A  table  Of  values  of  X,  X®,  Y  and  XY  is  formed: 


The  slope  of  the  least  squares  straight  line  is  determined  from  the  equation: 


D.  M>Anal2Sisi^fiJ|A^<<SMgle<iData 


1.  Coefficients  of  Fitted  Equations 


Sample 

EH 

mm 

B33I 

wvsm 

IX1X2Y 

mm 

Baal 

1 

-22.03 

-11.31 

70.53 

75.91 

-0.35 

112.79 

75.19 

2 

-18.68 

-13.92 

71.96 

-  72.86 

“4.79 

114.32 

76.21 

3 

-20.68 

-n.89 

68.96 

73.05 

-1.63 

112.23 

74.82 

4 

-16.76 

-  7.54 

72.88 

77.16 

-1.05 

117.05 

78.03 

5 

-I9.60 

-10.77 

72.42 

73.25 

-2.13 

112.26 

74.84 

6 

-21.22 

-10.26 

70.22 

73.08 

-1.1? 

114.11 

76.07 

7 

-21.76 

-  9.79 

75.54 

76.65 

-4.74 

115.69 

77.13 

8 

-20.46 

-10.53 

69.84 

73.47 

-3.26 

I09.96 

73.31 

9 

-20.8? 

-12.10 

64.51 

70.74 

-2.52 

106.48 

70.99 

10 

-22.55 

-  7.05 

70.01 

71.63 

-4.12 

UO.07 

73.38 

1 

m 

m 

K0| 

WB&m 

mu 

IMUl 

_  6 

MUMR 

■a 

d 

2 

1 

-3.67 

“1.89 

-0.09 

“2.33 

0.36 

— 

2 

-3.11 

-2.32 

-1.20 

-2.12 

-1.68 

3 

“3-45 

-1.98 

-o.4i 

-2.93 

-0.89 

15.02 

4 

-2.79 

-1.26 

-0.26 

-2.58 

-0.44 

15.02 

5 

-3.27 

-1.80 

-0.53 

-1.21 

-0.80 

13.81 

6 

-3.54 

-1.71 

-0.29 

-2.93 

-1.50 

15.63 

7 

-3.63 

-1.63 

-1.19 

-0.80 

-0.24 

13.55 

8 

-3.41 

-1.76 

-0.82 

-1.74 

0.08 

13.33 

9 

-3.48 

-2.02 

-0.63 

-3.24 

-0.13 

14.08 

L0 

-3.7 6 

-1.18 

-1.03 

-1.69 

-0.88 

13.94 

FittedJSguatioM 


1.  Y  =  13.85  -  3.67X1  -  1.89X2  -  2.33X12  +  0,36x/  -  o.09X1x2 

2.  Y  =  15.24  *  3.HXX  -  2.32X2  -  2.12Xi2 3 4 5 6 7 8 * 10  -  1.68x/  -  l,2Q3CjXg 

3.  Y  =  15.02  -  3.45Xi  -  1.98X3  -  2.93X12  -  O.89X22  -  0.41XjX2 

4.  Y  =  15.02  -  2.79Xi  -  1.26Xg  -  2.58x/  -  0.44.x/  -  0.26Xi3^ 

5.  Y  -  13.81  *  3.27X1  -  1.80Xg  -  I.21X12  *  O.80X22  -  0.53X1X2 

6.  Y  =  15.63  -  3.54Xi  -  I.7IX2  -  2.93X/  -  1.50X/  -  0.29XiX2 

7.  Y  -  13.55  -  3.63Xi  -  I.63X2  -  0.80X/  -  0.24X22  -  l.WXiXg 

8.  Y  =  13.33  -  3.4lXi  -  I.76X2  -  1.74X/  +  O.O8X22  -  0.823^X2 

9-  Y  -  14.08  -  3.48Xi  -  2.02Xg  -  3.24x/  -  0.13Xg2  -  0.63X^Xg 

10.  Y  =  13.94  -  3.76X1  -  1.18X2  -  I.69X/  -  O.88X22  -  i.OS^Xg 


com  cn  in  men  tn  co  cry 


3 .  .  Predicted  Points  froc  Fitted  Equations 


1)  13.75 

6)  13.32 

1)  12.32 

6)  12.42 

1)  6.23 

6) 

5.66 

2)  13.^3 

7)  15.70 

2)  11.24 

7)  11.68 

2)  4.81 

7) 

6.06 

3)  13.08 

8)  14.14 

3)  12.15 

8)  11.65 

_ 

3)  5.36 

8) 

5.68 

4)  13.79 

9)  12.80 

■ 

4)  13.32 

9)  11.93 

4)  7.69 

9) 

4.56 

5)  13.80 

10)  14.98 

5)  11.21 

10)  11.88 

5)  6.20 

10) 

5.4o 

1)  15.19 

6)  16.24 

1)  13.85 

6)  15.63 

1)  7.85 

6) 

9.16 

2)  16.23 

7)  16.38 

2)  15.24 

7)  13.55 

2)  10.01 

7) 

9.12 

3)  15.54 

8)  15.00 

- 

3)  15.02 

8)  13-33 

— 

3)  8.64 

8) 

8.18 

4)  15.23 

9)  14.32 

4)  15.02 

9)  14.08 

4)  9-65 

9) 

7.36 

5)  15.87 

10)  16.01 

5)  13.81 

10)  13.94 

5)  9.33 

10) 

8.49 

i)  17.35 

6)  16.16 

1)  16.10 

6)  15.84 

1)  10.19 

6) 

9.66 

2)  15.67 

7)  16.58 

2)  15.88 

7)  14.94 

2)  11.85 

7) 

11.70 

3)  16.22 

8)  16.02 

— 

3)  16.11 

8)  15.17 

— 

3)  10.14 

8)  10.84 

4)  15.79 

9)  15.58 

4)  15.84 

9)  15.97 

4)  10.73 

9) 

9.88 

5)  16.34 

10)  15.26 

5)  14.81 

10)  14.24 

5)  10.86 

10) 

9-82 

© 


© 


4.  Standard  Errors 


Bb 

rental 

10.19 

15.19 

mmm 

mem 

13.75 

nm  tti  1  ^ 

y 

16.94 

9.91 

15.77 

13.40 

7.71 

13.77 

11.92  6.6l 

TTyT 

6.m 

■»iBa 

0.5» 

0.14 

0.02 

0.40  0,38 

2 

1.8365 

V°*  54550 

-  0.739 

T 

--  U-5 

ITT 

2)(Y-y) 

O.56 

0.55 

0.43 

0.19 

0.60 

0.45 

O.39  0.07 

-i/i 

=  . 

Sggjlgj 

=•  0.72 

JO 

0.25 


0.35 

ysss 

3  ■ 


1.29  0.71  0.60  0.71 

-  ^1.16077  =  1.077 


0.45 


0.25 


3)(Y-y) 


1*.  StendarfiError^^ontinued^ 


4) 

(Y-y) 

0.66 

0.32 

0.3 

10  0.58 

0.29  0.02 

0.07 

0.04 

s* 

OKI 

■ 

-35180 

»  0.593 

5) 

(Y-y) 

0.65 

0.76 

O.s 

il  0.69 

1.20  0.40 

0.03 

0.46 

s= 

* 

Wl. 18293 

*  1.088 

^  3 

6) 

(Y-y) 

0.69 

0.45 

o.c 

>8  O.63 

0.55  0.15 

0.06 

0.08 

s= 

ss 

VO  *47230 

il  0.69 

«  0.687 

7) 

(Y-y) 

1.13 

>  3 

0.43 

o.e 

0.07  0.08 

0.46 

0.37 

s* 

^ 2.6698 

s 

VO -88993 

=  0.943 

8) 

(Y-y) 

0.77 

0.88 

0.88  0.28 

1.18  0.75 

0.46 

0.29 

s= 

^ 4.4707 

S 

V 1.49023 

=  1.221 

9) 

(Y-y) 

0.06 

0.23 

0.12  0.57 

O.67  0.05 

0.50 

0.44 

s= 

VH8* 

= 

Vo.  43027 

«  O.658 

10) 

(Y-y) 

0.13 

0.39 

1.39  0.95 

0.44  0.88 

0.82 

0.03 

s= 

fiH 

= 

V 1.54830 

p  1.244 
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5.  Comparison  of  Standard  Deviations  end  Standard  Errors 


■ 

mmm 

Variances- 

Variance 

Ratio 

Hi 

n 

1 

0.739 

1.054 

0.54550 

1.11022 

Ho 

0.92198 

0.960 

2 

0.720 

0.995 

0.51887 

0.98978 

No 

0.83281 

0.913 

3 

1.077 

0.883 

1.16077 

0.69433 

No 

0.84981 

0.922 

4 

0.593 

1.085 

0.3518c 

1.17772 

0.30 

No 

0.90241 

0.950 

5 

1.088 

1.129 

1.18293 

1.27372 

0.93 

No 

1.24346 

1.115 

6 

0.687 

0.788 

0. 47230 

0.62088 

0.76 

No 

0.57135 

0.756 

7 

0.943 

1.488 

0.88993 

2.21483 

0.40 

No 

1.77320 

1.332 

8 

1.221 

0.546 

1.49023 

0.29798 

5.00 

No 

0.69540 

0.834 

9 

0.656 

1.012 

0.43027 

1.02458 

0.42 

No 

0.82648 

0.909 

10 

1.244 

0.739 

1.54830 

0.54662 

2.83 

No 

0.88051 

0.938 

*  F  Ratio  Test  for  Comparison  of  Variance 


Reject  the  hypothesis  that  the  variances  are  equal,  at  the  95#  level 
of  significance,  if  the  variance  ratio  falls  outside  the  region  of  6.60  to 
0.068  /  \  Ho  means  no 


C.  ft  f  \  =  6*60,  ~  0 .068 

•975(3,6  d.f.)  F. 975(6,3  d.f.) 


significant  difference,  permitting  pooling  of  the  variances. 


**  Calculations  for  Pooling  Variance,  Sample  1 

Sp2  =  3(a.f.)  x  0.5^550  +  6(d.f .)  x  1.11022  «=  0. 92198 

9  (d.f.) 

Sp  =V  0.92198  =  0.960 
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6.^Reli^llit^^omifier^<iE2uation^ 

'  Based  on  Y-KS  where  K  =  5*4l4  for  N  =  15-6  *9  and  for  reliability J 

standards  of  99  «S %  at  95$  confidence .  J 

1.  3.67X2  +  1.89X2  +  2.33X22  -  O.36X22  +  0.09X2X2  =  13.85  -  5.20  -10.00—1.35 


2. 

3.11Xi 

+  2.32X2  +  2.12Xi^ 

+  l.68Xg2  +  1.20X1X2  =  15.24  -  4.94  -10.00=  0.30 

3. 

3. ^5Xi 

+  1.98X2  +  2.93Xi2 

+  O.89X22  +  O^^Xg  -  15.02  -  4.99  -10.00=  0.03 

4. 

2.79Xi 

+  I.26X2  +  2.58X12 

+  o.llXg2  +  o.26x,Xg  =  15.02  -  5.14  -10.00—0.12 

5. 

3.27X2 

+  I.8OX2  +  1.2IXi2 

+  O.8OX22  +  O.53X  Xg  =  13.81  -  6.04  -10.00—2.23 

6. 

3.5^ 

+  I.7IX2  +  2.93X12 

+  1.50X22  +  0.29X  Xg  =  15.63  -  4.09  -10.00=  1.54 

7« 

3.63X2 

+  1.63X£  +  O.SoX^2 

+  0.24x22  +  I.19X  Xg  =  13.55  -  7.21  -10.00—3.66 

2  2 

8.  3.1HX2  +  1.76X2  +  1.74X  -  0.08x2  +  0.82XiX2  =  13.33  -  4.52  -10.00— 1. 19 

9.  3A8Xl  +  2.02X2  +  3.24X  2  +  0.13X22  +  O.63XX  =  14.08  -  4.92  -10.00—0.84 

10.  3.76X,  +  1.18X  +  1.69X  2  +  O.88X  2  +  1.03X  X  =  13.94  -  5.08  -10.00— 1.14 

1  2  1  2  1  2 
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E.  _Anal2si8_o£j]B^iSamglei>DateM_ 

1.  Coefficients  of  Fitted  Equations  for  "B"  Samples 


Sample 

SXiY 

XXgY 

@3 

^Y 

2/32Y 

1 

-40.74 

-19.91 

141.40 

150.03 

-0.93 

225.00 

150.00 

2 

-4o.4l 

-24.62 

141.73 

143.54 

-8.13 

226.21 

150.81 

3 

-40.41 

-19.07 

143.81 

147.79 

-5.43 

229.77 

153.18 

4 

-35.75 

-18.67 

143-31 

151.81 

-4.83 

231.04 

154.03 

5 

-41.91 

-20.29 

140.03 

142.61 

218.29 

145.53 

6 

-38.88 

-15.99 

142.86 

145.63 

-0.97 

228.41 

152.27 

7 

-43.87 

-19.89 

147.07 

151.27 

-7.64  1 

230.66 

153.77 

8 

-42.52 

-22.98 

141.68 

146.74 

-6.19 

220.39 

146.93 

9 

-38.35 

-22.80 

140.35 

146.90 

-8.56 

224.63 

149.75 

10 

-41.68 

-16.28 

140.40 

145.04 

-7.13  _ 

223.66 

149.11 

| 

m 

m 

HU 

WHOM 

MBIM 

12 

12 

'B:  

5 - 

5 

1 

-3.40 

-1.66 

-0.12 

-2.15 

0.01 

13.93 

2 

-3.37 

-2.05 

-1.02 

-2.27 

-1.82 

15.29 

3 

-3.37 

-1.59 

-0.68 

-2.34 

-1.35 

15.23 

4 

-2.98 

-1.56 

-0.60 

-2.68 

-0.56 

14.99 

5 

-3.49 

-1.69 

-0.58 

-1.38 

-0.73 

13.53 

6 

-3.24  , 

-1.33 

-0.12 

-2.35 

-1.66 

15.37 

7 

-3.66 

-1.66 

-0*96 

-1.68 

-0.63 

14.35 

8 

-3.54 

-1.92 

-0.77 

-1.31 

-0.05 

13.15 

9 

-3.20 

-1.90 

-1.07 

-2.35 

-0.71 

14.52 

10 

-3.47 

-1.36 

-0.89 

-2.16 

-1.02 

14.56 

Fittedjg^tions 


1*  Y  -  13.93 

2.  Y  =  19.29 

3.  Y  =  15.23 

4.  Y  =  14.99 

5.  Y  =  13.53 

6.  Y  =  15.37 

7.  Y  =  14.35 

8.  Y  -  13.15 

9.  Y  =  14.52 
10.  Y  =  14.56 


3.4ox1  -  1.66X2 
3.37Xx  *  2.05X3 
3.37XX  -  1.59X2 
2.98X1  *  1.56X2 
3.49X!  -  I.69X2 
3.24XX  -  I.33X2 
3.66X1  *  I.66X2 
3.54Xi  -  I.92X2 
3.20Xi  *  1.90X2 
3.47X1  -  1.36X2 


2.15X12  +  O.OIXg2 
2.27Xi2  -I.82X22 
2.34Xi2  -  I.35X22 
2.68Xi2  *  O.56X22 
1.3SXi2  -  O.73X22 
2.35X12  -  1.66Xg2 
I.68X12  -  O.63X22 
1.31XJ2  -  0.05Xg2 
2.35X12  -  O.TLXg2 
2.l8Xi2  -  1.02Xg2 


0.^X2 

I.02X1X2 

O.68X1X2 

0.60X1X2 


0.89X1X2 
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12.31 

14.15 

6.81 

9.33 

6.75 

4.38 

5-36 

8.48 


^Linear  Terms:  (2X3Y) 

12 

♦^Quadratic  Terms:  (fX^Y 


+  (sXgSy  -  2/35Y)2  +  (SX^XgY )‘ 


***Are  quadratic  terms  significant  at  the  95#  level  of  significance? 

To  he  significant  the  ratio  must  exceed.  F#^  g  12  *  3*^9 

the  lack  of  fit  term  significant,  at  the  95#  level  of  significance, 
in  comparison  to  the  experimental  error?  To  he  significant  ®  h  0  rrftio 
must  exceed  j^q  *  3.64.  §3  4,9 


3.  Predicted  Points  from  "B"  Sample  Pitted 


tions 


1)  13.65 

2)  13.0U 

3)  lb.00 
b)  13.77 
5)  13.80 


11.11 


6) 

6.67 

7) 

0.76 

8) 

0.06 
ef  06 

y ) 
10) 

0.6b 

1)  10.18 

2)  16.39 

3)  16.26 
b)  10.20 
0)  10.6U 


2)  10.29 
3)  10.23 


6) 

10.37 

7) 

lb.  30 

8) 

13.10 

9) 

lb.02 

10) 

lb.  06 

8.38 

9.60 

9.02 

9.33 


6) 

9.78 

7) 

9.01 

8) 

8.30 

9) 

8.97 

10) 

8.91 

1) 

16.73 

2) 

10.60 

3) 

10.82 

b) 

10.69 

0) 

16.02 

1)  10.17 

2)  10.00 

3)  lo.llb 

.93 

0)  10.20 


6)  9.07 

7)  11.00 

8)  10.9b 

9)  11.23 
10)  10.1b 


le  Standard  Errors 


>.60  10.17  10.18  13.93  8.38  13.60  12.28  6.61 


Ai-y) 

N-0 

■Vs 

.6208 

11 

*  yjo.lmo  7 

-  0.6U8 

- 

0.92 

O.bO 

0.09 

0.2b  0.96 

0.3b 

0.07 

0.02 

1.16 

0.08 

O.b0 

1.20  1.82 

0.10 

0.18 

0.60 

3) 

(i-y) 

0.8? 

o.o5 

0.57 

0.92 

1.1*8 

1.63 

0.31 

0.29 

1.13 

0.50 

0.53 

0.1*1* 

0.82 

0.82 

0.39 

0.38 

s.  y 

/10.63SL 
'  U 

■1 

VO.  96692 

•  O.983 

U) 

(i-y) 

0.81 

0.12 

0.36 

0.61 

0.03 

0.00 

0.52 

1.01* 

0.27 

0.16 

0.71* 

1.1*0 

0.39 

0.63 

0.81* 

1.70 

B 

rrer 

31 

•  ^/o.8U325 

-  0.918 

5) 

(i-y) 

0.33 

1.1*2 

0.7I* 

0.97 

0.53 

0.1*0 

0.07 

1.00 

0.5k 

0.1*8 

0.35 

0.28 

1.09 

0.10 

0.1*5 

0.30 

S  *  ^ 

/  7.2731 
11 

-  -y/o.66119 

*  0.813 

6) 

(i-y) 

1.1*9 

0.36 

0.06 

0.37 

0.07 

0.08 

0.02 

1.09 

1.22 

0.58 

0.99 

0.98 

2.31 

0.16 

1.07 

0.26 

ra 

• 

■t 

V  1.2751a 

’  -  1.129 

7) 

(i-y) 

0.69 

0.27 

0.66 

1.1*9 

0.18 

0.78 

0.81* 

0.67 

1.17 

1.13 

0.5U 

1.1*9 

0.29 

0.09 

1.01 

0.36 

n 

1  II.399U 

U 

• 

Vl.03631 

-  1.018 

8) 

(i-y) 

0.62 

0.78 

0.50 

0.1*6 

1.30 

0.79 

0.93 

o.la 

1.26 

0.66 

0.1*0 

0.1*6 

0.39 

0.69 

1.56 

0.36 

S'\ 

/ 10.3877 

1  n 

■ 

Vo.9l*l*3l* 

.  0.972 

9) 

(I-y) 

0.32 

1.58 

1.17 

1.01 

0.91* 

1.08 

0.52 

1.15 

1.15 

1.06 

0.89 

1.37 

2.65 

1.77 

0.Q5 

0.02 

D 

/23.90W. 

31 

t 

P 

^2.17310 

’  -  1.1*71* 

*0) 

(I-y) 

0.79 

0.71 

1.23 

0.33 

0.02 

1.50 

1.12 

0.27 

0.80 

0.32 

0.82 

0.28 

1.96 

0.1*3 

1.01* 

0.70 

HI 

p 

■ 

•^1.21991 

‘  -  iaoU 
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5.  Comparison  of  "B"  Sample  Standard  Deviations  and  Standard  Errors 


Sample 

Std. 

Error 

Std. 

Dev. 

—  Variance  — 

Variance 

Ratio 

I 

^9 

SP 

(Std.Error)2 

(Std. Dev. )2 

O.6U8 

0.887 

0.U2007 

0.786U8 

0.53 

No 

0.6185k 

0.786 

0.935 

0*989 

0.8?kk5 

0.97779 

0.89 

No 

0.930k3 

0.965 

0.983 

0.665 

0.96692 

0.kk202 

2.19 

No 

0.68260 

0.826 

0.918 

1.021 

0.8U325 

1.0U261 

0.81 

No 

0.9512k 

0.975 

5 

0.813 

1.067 

0*66119 

1.13809 

0.58 

No 

0.91951 

0.959 

6 

1.129 

1.108 

1.275U1 

1.22775 

1.0k 

No 

1.2k959 

1.118 

7 

1.018 

1.313 

1.03631 

1.723U2 

0*60 

No 

I.k08k9 

1.187 

8 

0.972 

0*830 

0.9kk3k 

0*68952 

1*37 

No 

0.80631 

0.898 

H 

1.U7U 

0.915 

0.8378k 

No 

l.kk983 

l*20k 

WM 

1.10k 

0.708 

2.UU 

No 

0.83035 

0.911 

*  F  Ratio  Test  for  Comparison  of  Variance 


Reject  the  hypothesis  that  the  variances  are  equal,  at  the  95?  level  of 
significance,  if  the  variance  ratio  falls  outside  the  region  of  3.20  to 


0.293 


(F 


3.20  , 


.975(11,13  d.f.)  -.97?(13>u  d.f.) 

no  significant  difference,  permitting  pooling  of  the  variances. 


0.293).  No  means 


6.  J^^Saagle^^Reliabili^^Boundar^^Eguations^ 

(Based  on  Y  -  KS  where  K  -  U.171  for  N  -  30  -  6  •  2k  and  for  relia-J 
bility  standards  of  99.9?  at  95?  confidence.  J 

1,  3.1^  +  1.66Z2  +  2.15X!2  -  O.OIXg2  +  O.lZXjXg  -  3.93  -  3.28  -  10.00-0.65 

2.  3.37X1  +  2.05X3  +  2.27X12  +  1.82X22  +  I.O2X3X2  -15.29  *  k,03  -  10.00-1.26 

3.  3.37X^  +  1,59X2  +  2,3hXj2  +  1.35Xg2  +  0.68X3X2  -15.23  -  3.U5  -  10.00*1,78 

k  2.98X1  +  1.56Xg  +  2.68X12  +  0.56Xg2  +  0.60X3X2  -ik.99  -  U.07  -  10.00-0.92 

5.  3.1*9X1  +  I.69X2  +  1.38X12  +  0.73Xg2  +  0.58X3X2  -13,53  -  k.00  -  10.00--0.U7 

6.  3.2UX^  +  1.33Xg  +  2.3pX^2  +  1.66Xg2  +  0.12X^Xg  -15.37  -  U.66  -  10.00-0.71 

7.  3.66X1  +  1.66X2  +  I.68X12  +  0.6^2  +  0.96X1X2  -1U.35  -  U.95  -  10.00-0.60 

8.  3.5UX1  +  1.92X2  +  1.31X32  +  O.ogXg2  +  0,77^Xg  -13,15  -  3.75  -  10.00-0.60 

9.  3.2CQ^  +  1.90X2  +  2.35X32  +  0.71X22  +  I.O7X1X2  -1U.52  -  5.02  -  10.00-0.50 

10.  3.U7Xi  +  1.36X2  +  2.183^2  +  1.02X22  +  0.89XiXg  -1U.56  -  3.80  -  10.00-0.76 

'^"Sample  Reliability  Boundary  Equation  for  Infinite  Number  a t  Sample 

3.21^  +  1.69X2  +  2.13X12  +  0.9UXg2  +  OJi^Xg  -  lk.61  -  k.17  -  tt.00  -  O.UU 
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is  or  C 


•  Coefficients  of  Fitted  Eouations 


Sample  I 


B 
8 
3 

-75.9k 
-81.5k 
-7k. 60 
-8k. 72 
-82.79 
-75.80 
-79.27 


Sample 


SXgY 

*Xi2Y 

mmmm 

SY 

2/3 

-kO.56 

287.81 

301.68 

-  2.k5 

k58.58 

305.72 

-kk.16 

285.  kk 

29k. 18 

-17.66 

k55.90 

303.93 

-39.76 

290.69 

297.82 

-  7.05 

k6o.27 

306.85 

-38.61 

288.26 

302. k7 

-  7.53 

k59.29 

306.19 

-39.81 

280.58 

287.91 

-  7.29 

442.80 

295.20 

-33.5k 

287.32 

293.30 

-  0.79 

k56.?6 

30k. 51 

-38.78 

286.46 

29k. 98 

-  7.92 

k55-05 

303.37 

-kk.82 

28k. 63 

290.62 

-  7.02 

kk5.53 

297.02 

-kl.27 

280.56 

29k. 83 

-12.02 

450.87 

300.58 

-35.1k 

28k. 19 

29k. 38 

-9.59 

452.29 

301.53 

SUES! 

-fX^Y 

2h 

2k 

j.6 

-3.3k 

-1.69 

-0.15 

-3.3k 

-1.84 

-1.10 

-3.18 

-1.66 

-0.44 

-3.16 

-1.6l 

-0.47 

-3.k0 

-1.66 

-0*46 

-3.11 

-l.ko 

-0.05 

-3.53 

-1.62 

-0.50 

-3.k5 

-1.87 

-0.44 

-3.16 

-1.72 

-0.75 

-3.30 

-1.46 

-O.60 

5X2%  - 


-2.2k 

-2.31 

-2.02 

-2.2k 

-1.83 

-2.15 

-2.11 

-1.55 

-2.50 

-2.17 


-0.5f 

-1.22 

-1.13 

-0.k7 

-0.91 

-l.ko 

-1.05 

-0.80 

-0.72 

-O.89 


k.57 
5.02 
4.88 
k.56 
lk.13 
15.06 
lk.  75 
13.9k 

lk.67 

lk.6l 


Fitted  E 


1.  Y  ■  lk.57  -  3.3^2.  -  1.69Xg  - 

2.  Y  -  15.02  -  3.3kXx  -  1.8kX2  - 

3.  Y  -  lk.88  -  3.18X1  -  I.66X2  - 

k.  Y  -  lk.56  -  3.1%  -  I.6IX2  - 

5.  Y  -  lk.13  -  3.kOXi  -  1.66X2  - 

6.  Y  «  15.06  -  3.HX1  -  l.kOXg  - 

7.  Y  «  lk.75  -  3.53X1  -  1.62X2  - 

8.  Y  -  13.9k  -  3.k5Xi  -  1.87X3  - 

9.  Y  ■  lk.67  -  3.I6X1  -  1.72X3  - 

10.  Y  -  lk.6l  -  3.3OX1  -  1.46Xg  - 


-  2.31Xi^ 

-  2.02Xi2 

-  2.2kXi£ 


I.83X12 
2.15Xi2 
2. l^2 
1.55Xi2 
2.503^2 
2.17Xi2 


0.51X2*  -  0.15Xi3^ 
1.223^2  .  i.lOXiXg 

1.13X22  -  O.klXi^ 
0.k73^2  „  o.k730L3C2 
0.913^2  -  0.k6XiXg 
l.kOXg2  -  0. 053^3^ 

1.053^2  .  0.503CjX2 
0.80X32  -  0.hkX1X2 
0.72Xg2  -  0.75X1X3 
O.89X22  -  O.6OX1X0 


3.  "C^Ssggle^tmidirdJErroi^^ ^craitinued^ 


2)(Y-y) 

0.80 

1.28 

0.11 

0.36 

0.21 

0.27 

1*09 

0.$8 

0.21  0.03 
0.80  1.52 
0.70  O.83 

0.05  6.61 

1.25 

1.53 

0.59 

0.21 

0.21 

0.70 

1.05 

0.20 

1.11 

0.72 

1.58 

0.11 

0.15 

0.17 

O.96 

0.13 

3)(Y-y) 

S‘V 

0.91 

19.8010 

27 

0.l6 

-  ^6*73318 

0.77  0.58 

=  0.856 

1*61 

1*31 

O.5O 

0.81 

1.08 

0.29 

6*73  0.78 

0.66 

1.U 

O.58 

0.17 

0.19 

0.25 

6.I3  0.12 

1*36 

1.52 

0.02 

1.12 

0.36 

0*10 

0*11  O.ll 

0.78 

1.69 

1.50 

O.98 

S'V 

0.52 

22.9175 

27 

0.29 

*  ^/o  *81991 

0.56  l.Ol 

-  0.922 

0.21 

0.11 

0.91 

1.05 

*>(i*y) 

0.56 

0.33 

0.51  0.97 

0*57 

0.71 

1.23 

1.69 

0.7! 

0.89 

1.39  O.ll 

l.lo 

1.17 

0.77 

0.28 

0.51 

0.57 

0.I3  0.66 

0.15 

0.80 

1.09 

0.65 

s.  ^ 

0.71 

21.1317  *  ^.79388  - 

27 

1.52  0.68  O.37 

0.891 

0.77 

0.18 

0.37 

0.80 

5)(Y-y) 

0.92 

O.38 

O.ll  0.32 

1.33 

0.32 

O.89 

0.50 

0.19 

O.56 

0.22  O.Ol 

0.53 

0.9I 

1.11 

O.83 

0.03 

0.39 

0.92  1.01 

0.95 

1.51 

0.6l 

O.26 

6)(Y-y) 

s  * 

1.1*3 

/IS.  9592 
1  "  27 

0.63 

*  ^f0. 70219 

0.31  0.03 

*  0.838 

0.08 

0.21 

0.11 

1.36 

1.23 

0.85 

0.7I  0.66 

2.30 

0.29 

1.20 

6.53 

0.62 

1.03 

1.06  0.95 

1.77 

0.60 

2.08 

0.62 

1.65 

0.65 

0.31  0.83 

1.76 

0.02 

2.02 

0.22 

s  » 

1 

/ 37.7136 
1  27 

-  V1*  39791 

=  1.183 
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3.  iJ[|Ci[|ii^emgle_iStandard^iErrors_£contimea21 


7)(Y-y) 

0.75 

1.09 

0.83 

1.88 

O.09 

1*78 

0.86 

0.49 

1.23 

0.31 

0.37 

1.10 

O.38 

1.09 

1.03 

0.18 

0.62 

1.08 

0.08 

0.37 

0.06 

0.23 

0.75 

1.82 

0.18 

0.64 

1.22 

1.72 

O.O5 

1.13 

0.18 

1.16 

s  - 

J28.6835 

Vi .06235 

-  1.031 

T  27 

0.61 

1.26 

0.03 

0.33 

1.95 

0.23 

0.84 

0.13 

8)(Y-y) 

1.27 

0.18 

0.07 

1.25 

0.26 

1.25 

1.65 

0.08 

0.37 

0.01 

1.09 

1.31 

0.47 

0.32 

0.37 

1.00 

0.28 

0.88 

0.56 

0.51 

0.16 

1.48 

0.33 

0.42 

s  = 

J22.6375 

b 

^0. 83843 

«  0.911 

«r 

D 

T  27 

9)(Y-y) 

0.36 

1.11 

1.13 

1.16 

O.98 

0.89 

0.20 

1.68 

1.11 

1.53 

0.85 

1.22 

2.6l 

1.96 

0.37 

0.55 

0.71 

0.06 

1.67 

0.97 

0.94 

0.44 

0.22 

0.22 

1.64 

0.11 

0.82 

1.21 

0.14 

0.84 

1.36 

O.83 

s  « 

J38.9895 

t? 

^1.44398 

*  1.202 

*  27 

10)(Y-y) 

1.07 

0.89 

1.12 

0.29 

0.21 

1.88 

1.18 

0.8l 

0.52 

0.50 

0.71 

0.32 

1.73 

0.05 

0.98 

1.24 

0.98 

0.69 

0.77 

0.26 

1.35 

0.78 

0.30 

0.72 

0.15 

0.60 

0.44 

0.96 

O.65 

0.18 

0.67 

0.31 

s  «= 

1 

B 

^0.857^6 

=  0.926 
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4.  Comparison  of  "C"  Semple  Standard  Deviations  and  Standard  Errors 


■ 

BH 

M 

Variance 

1 

HI 

■ 

KSShH 

O.805 

0.715 

0.64799 

0.51193 

1.27 

Ho 

0.57996 

0.762 

0.356 

0.834 

0.73348 

0.69610 

1.05 

Ho 

0.71479 

0.845 

0.922 

1.043 

0.84991 

1.08686 

0.78 

No 

0.96839 

0.984 

0.891 

0.921 

0.79388 

0.84844 

0.94 

Ho 

0.906 

0.838 

0.985 

0.70219 

0.97113 

0.72 

Ho 

0.83666 

0.915 

6 

1.183 

1.035 

1.39791 

1.07163 

1.30 

HO 

1.23477 

1.111 

B 

1.031 

1.164 

1.06235 

1.35407 

0.78 

No 

1.20821 

1.099 

K9 

0.916 

0.850 

0.83843 

0.72324 

1.16 

No 

0.78109 

0.884 

KB 

1.202 

0.917 

1.44398 

0.84103 

1.72 

No 

1.14251 

I.069 

ESI 

0.926 

0.866 

0.85766 

0.75011 

1.14 

No 

0.00369 

0.897 

*  F  Ratio  Test  for  Comparison  of  Variance 


Reject  the  hypothesis  that  the  variances  are  equal,  at  the  95$  level 
of  significance,  if  the  variance  ratio  falls  outside  the  region  of  2.17  to 
0.461 

(*F.975(27,27  a.f.)  +  2*1^  — _i - - -  ■  0.461)1 

V  .975(27,27  d.f.)  ' 

Ho  means  no  significant  difference,  permitting  pooling  of  the  variances. 
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5.  "C"  S eerie  Reliability  Boundary  Equations 

Based  on  Y-KS  where  K*  3  •731  for  N*  60  -  6  * 
standards  of  99*9#  at  95#  confidence. 


2  2 

1.  3.34x,  +  1.69X  +  2.24x  +  0.51X  +  0.15X  x  * 

1  2  l  2  12 

2  2 

2.  3.3IX  +  1.84X  +  2.31X  +  1.22X  +  1.10X  X  = 

1  2  1  2  12 

2  2 

3.  3.18X  +  1.66X  +  2.02X  +  1.13X  +  0.44X  X  * 

1  2  1  2  1  2 

2  2 

4.  3.16X  +  1.61X  +  2.24X  +  0.17X  +  O.lTX  X  * 

1  2  1  2  12 

2  2 

5.  3»*»OX  +  1.66X  +  I.83X  +  0.9IX  +  0.46x  X  = 

1  2  1  2  12 

6.  3.HX  +  1.40X  +  2.15X  2  +  1.40X  2  +  0.05X  X  * 

1  2  1  2  12 

7.  3.53X  +  1.62X  +  2.11X  2  +  1.05X  2  +  0.50X  X  = 

1  2  1  2  1  2 

8.  3.45X  +  1.8?X  +  I.55X  2  +  0.80X  2  +  0.44x  X  = 

1  2  !  2  12 

9.  3.16X  +  1.72X  +  2.50X  2  +  0.72X  2  +  0.75X  X  * 

1  2  1  2  12 

2  2 

10.  3.30X  +  1.46X  +  2.17X  +  0.89X  +  0.60X  X  * 

1  2  1  2  12 


if 


54  and  for  reliability^ 

14.57  -  2.84  *10.00*  1.73 
15.02  *  3.15  -10.00*  1.87 
14.88  *  3.67  -10.00*  1.21 
14.56  *  3.38  -10.00*  1.18 
14.13  *  3.41  *10.00*  0.72 
15.06  *  4.15  *10.00*  0.91 
14.75  *  4.10  *10.00*  O.65 
13.94  *  3.30  -10.00*  0.64 
14.67  -  3.39  -10.00-  0.68 
l4.6l  -  3.35  -10.00*  1.26 


"C"  Sample  Reliability  Boundary  Equation  for  Infinite  Humber  of  Samples 

3.24X  +  1.69X  +  2.13X:  2  +  0.94x  2  +  O.43X  X  *  l4.6l  -  3.73  -10.00*0.88 
1  2  1  2  12 
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